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The  process  developed  produces  double  doped,  both  gallium  and 
boron  substrates. 

The  characteristics  of  the  gallium  Impurity  center  in  silicon  have 
been  determined  by  measurements  on  the  MOSFET  device  structure,  and 
specifically  the  thermal  emission  rate  of  holes  from  neutral  gallium 
centers  determined.  To  the  best  of  the  author's  knowledge  the  emission 
rate  of  gallium  has  not  been  previously  measured.  A very  significant 
field  enhancement  has  been  observed  in  this  thermal  emission  rate  and 
is  described  by  the  Poole-Frenkel  effect. 

Measurements  have  been  made  of  the  characteristics  of  the  gallium 
doped  MOSFET  (IRFET)  as  an  infrared  detector.  It  has  been  found  that  the 
device  operated  according  to  previously  published  models  and  design 
criteria  and  in  a manner  consistent  with  that  previously  observed  for 
the  indium  and  gold  doped  devices. 

These  results  demonstrate  that  the  gallium  doped  infrared  sensing 
MOSFET  (IRFET)  will  operate  in  the  far  infrared,  8 to  14  micrometer 
wavelength  range,  and  as  such  might  be  particularly  suited  to  large 
scale  integrated  infrared  imaging  arrays  and  applications. 
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TECHNICAL  REPORT  SUMMARY** 

The  gallium  diffusion  process  from  Emulsitone  doped  oxide  sources 
was  developed  using  MOS  capacitor  structures.  The  process  was  then  modi- 
fied and  applied  to  MOSFET  device  structures  by  doing  the  gallium  diffu- 
sion after  or  during  gate  oxidation. 

These  latter  MOSFET  devices  were  then  characterized  in  a low  tempera- 
ture cryo-tip  in  the  range  20  to  50  °K.  A rather  significant  field  en- 
hancement was  observed  in  the  thermal  emission  rate  of  the  gallium  center 
in  silicon,  the  longest  time  constant  observed  being  of  the  order  1.0 
second  at  20  °K.  Detailed  measurements  were  made  of  the  change  in  tran- 
sistor characteristics  due  to  thermal  emission  or  optical  emission  of  holes 
from  the  gallium  centers  due  to  incident  infrared  radiation.  The  response 
characteristics  correspond  to  the  model  previously  developed  for  the  infra- 
red sensing  MOSFET(IRFET)  and  are  consistent  with  previous  measurements  on 
indium  and  gold  doped  devices. 

Responsivities  of  the  order  of  25  milliwatts/microjoule  have  been 
readily  observed  in  non-optimum  device  structures  and  the  transistor 
characteristics  have  been  found  to  be  relatively  insensitive  to  temperature 
in  the  range  20  to  30  °K. 

More  work  remains  to  be  done  in  obtaining  more  optimum  device  structures, 
detailed  spectral  response  characteristics,  and  investigating  uniformity  of 
response. 


**  This  report  is  based  on  an  MSEE  Thesis  submitted  to  the  Department  of 
Electrical  Engineering,  University  of  Arkansas,  by  Hamman  ElAbd. 
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I.  INTRODUCTION 


This  thesis  describes  the  effect  of  inf rnred  radiation  on  the 
gallium  impurity  centers  in  the  surface  depletion  region  of  a metal-oxide- 
silicon  field  effect  transistor  (MOSFET). 

The  research  objectives  were  to  measure  the  thermal  emission  rate 
dependance  on  temperature,  the  charge  state  device  characteristics  and 
the  optical  emission  rate-wave  length  dependance  in  the  far  infrared 
region  (8  to  14  micrometer  wavelength). 

Gallium-doped  MOSFETs  were  designed  and  fabricated  in  the  University 
of  Arkansas  Solid-State  Devices  Laboratory.  The  fabrication  processes 
for  the  gallium  doped  MOS  capacitor  and  MOS  field  effect  transistor  are 
described  in  the  first  chapter  of  this  thesis.  We  did  not  spend  much  time 
on  the  Ga.  doped  MOS-C,  but  used  it  mostly  to  determine  if  we  were  getting 
enough  gallium  into  the  fabricated  device,  i.e.  the  MOS-C  was  used  as  a 
process  monitor.  In  other  words,  the  objective  was  to  compare  the  theo- 
retical and  experimental  characteristic  of  the  high  frequency  capacitance- 
voltage  (C-V)  curves,  and  deduce  from  this  the  total  impurity  concentra- 
tion at  the  surface  of  the  wafers.  The  gallium  and  boron  concentrations 
were  determined  from  the  low  temperature  measurements  on  the  gallium- 
boron  doped  MOSFETs. 

The  key  issue  in  the  fabrication  process  of  Ga-B  doped  MOSFET  was 
to  ensure  satisfactory  gallium  concentration  in  the  end-device  by  pre- 
venting the  gallium  atoms  from  escaping  from  the  Si-surface  into  the 
atmosphere  during  predeposition  or  thermal  oxidation. 

The  thermal  and  optical  response  of  the  gallium  centers  could  only 
be  observed  at  low  temperatures.  Satisfactory  time  constants  were  measured 
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in  the  temperature  range  21°K  to  50°K.  The  device  was  mounted  in  a 
liquid  hydrogen  cooled  refrigerator.  Operation  at  the  lower  end  of  the 
temperature  range  was  only  possible  when  the  thermal  load  on  the  cryotip 
was  reduced  to  a minimum. 

To  minimize  the  300°K  background  radiation  which  peaks  in  the  device 
response  range,  the  IR  filter  was  placed  inside  the  cryogenic  refrigerator 
(cold  filter  operation).  The  IR  filter  was  attached  to  the  radiation 
shield  of  the  cryogenic  refrigerator,  which  is  held  at  liquid  nitrogen 
temperature. 

The  thermal  and  charge  state  device  characteristics  enabled  us  to 
preliminarily  study  the  electric  field  effect  of  the  re-emission  of  the 
trapped  holes  in  the  gallium  centers.  This  effect  is  well  explained  by 
a lowering  of  the  Coulomb  potential  barrier  about  the  gallium  centers 
18, 9)  (the  Poole-Frenkel  effect). 

The  operation  of  the  infrared  sensing  MOSFET  (IRFET)  can  be  analyzed 
using  the  well  known  MOSFET  equations  [2].  The  experimental  results  were 
found  in  most  cases,  to  verify  the  model  equations  and  with  an  accuracy 
sufficient  for  most  engineering  applications.  Although  the  device  is  de- 
signed to  operate  as  a photon  detector  in  the  8 to  14  micrometer  wavelength 
range,  the  activation  energy  of  the  gallium  center;  calculated  from  the 
thermal  response  data,  suggests  that  the  gallium  center  should  respond 
to  photon  energies  from  1 micrometer  in  the  near  infrared  and  up  to  41 
micrometer  wavelength  in  the  far  infrared  range. 

The  double  doped  Si-IR  sensing  MOSFET  is  characterized  by  negative 
photoconductivity,  high  responsivlty , Integrating  function  and  static 
ROM  capability. 


TABLE  I ■ INFRARED  DETECTORS, 

TRANSMISSION  MATERIALS  AND  SOURCES 

ELECTRONICS  February  1972 


Considerable  amount  of  research  has  recently  been  reported  on  ex- 
trinsic silicon  infrared  detectors  [ 1 , 2, 3,4 , 5, 6, 7 ] because  of  the 
obvious  benefits  obtained  by  employing  the  silicon  technology.  Table 
shows  tha  available  IR  detectors,  transmission  materials  and  sources 
and  their  location  in  the  IR  spectrum. 
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II.  GALLIUM  DOPING  BY  DIFFUSION  FROM 
DOPED  OXIDE  SOURCES 

Gallium  diffusion  differs  from  the  conventional  diffusion  process 
(such  as  boron  and  phosphorus)  encountered  in  semiconductor  manufac Luring 
in  that  gallium  out-dif fuses  rapidly  when  silicon  is  thermally  oxidized, 
and  silicon  dioxide  does  not  mask  gallium  [10]  . The  diffusion  constant 
of  gallium  in  Si02  is  estimated  to  be  a factor  of  625  greater  than  its 
diffusion  constant  in  silicon  at  1200°C  [10]  . The  Si02  is  therefore 
transparent  to  the  gallium  for  practical  purposes.  The  diffusion  con- 
stant of  gallium  in  silicon  is  approximately  four  times  that  of  boron 
at  the  same  temperature. 

The  objective  of  this  section  of  the  study  was  to  establish  a process 
for  producing  double  doped  (gallium  and  boron)  n-channel  MOSFET’s.  The 
approach  was  to  begin  with  p-type  borop  wafers  and  diffuse  phosphorus  for 
the  source  and  drain  regions  and  gallium  in  the  channel  region  under  the 
gate.  Gallium  doping  by  diffusion  from  doped  oxide  sources  was  carried 
out  at  the  last  phase  of  the  MOSFET  fabrication  process.  A liquid  diffusion 
source  (gallium  silicafilm)  is  spun  over  the  entire  wafer.  The  wafer  is 
allowed  to  dry  then  subjected  to  a predeposition  drive-in  cycle.  It 
was  found  that  the  gallium  is  violatile  at  high  temperatures  and  the 
wafers  must  be  stacked  during  predeposition  to  prevent  out  diffusion. 

According  to  the  manufacturer  [ll],  the  source  produces  a co-error 
function  distribution  for  predeposition  time  up  to  one  hour.  After 
predeposition,  the  gallium  silicafilm  is  removed  by  etching  and  a drive-in 
step  is  performed  to  bring  down  the  surface  concentration  somewhat 
and  produce  a deeper  more  uniform  diffusion  profile. 
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To  insure  that  gallium  does  not  invert  the  source  and  drain  regions, 

a phosphorus  source  with  a surface  concentration  higher  than  the  solid 

solubility  of  gallium  is  used.  The  phosphorus  source  concentration  was 
21  3 

10  /cm  , which  is  higher  than  the  maximum  value  of  gallium  solid 

solubility  A x lO^cm-^  observed  at  1250°C.  Two  different  gallium 

1 9 -3  17  a 

source  concentrations  were  used  10  cm  and  10  cm-  . The  resulting 
total  doping  concentration  was  2.25  x lO^cm--*.  The  calculated  gallium 
concentration  was  0.A5  x lO^cm-^. 

The  MOSFET's  were  produced  on  2 ohm-cm  (100)  orientation  boron 
doped  substrate.  The  channel  W/L  ratio  is  8.17.  The  gate  oxide  estimated 
by  viewing  the  oxide  perpendicularly  under  white  light  and  varified  by 
the  C-V  measurements  was  around  7000A°. 

The  section  below  describes  the  experiments  performed  and  the  results 
in  establishing  the  gallium  process. 

A.  VERIFICATION  OF  GALLIUM  DOPINC 

Many  MOS  capacitors  were  fabricated  on  P and  N-type  silicon  wafers. 
Gallium  was  predeposited  on  the  wafer  surface  for  various  periods  of 
time.  Also  the  drive  in  diffusion  period  was  varied.  The  oxide 
thickness  was  varied  between  2000A°,  3000A°,  and  4000A°.  The  gallium 
diffusion  in  the  Si  wafer  was  observed  with  satisfactory  concentrations. 

The  object  wafer  was  stacked  between  two  buffer  n-type  wafers 
during  the  predeposition  and  drive-in  periods,  to  prevent  gallium  tiom 
escaping  the  surface  during  heat  treatment. 

IQ 

Gallium  diffusion  was  performed  from  a 10  cm-J  source  into  a 
group  of  10^cm“3  phosphorus  doped  N-type  Si  wafers.  The  result  of  this 
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experiment  was  a definite  surface  inversion;  verified  by  measuring  the 
capacitance  characteristics  of  the  capacitor  versus  voltage.  The  C-V 
characteristics  corresponded  to  p-type  surface  instead  of  an  N-type 
surface,  see  Fig.  1. 

From  the  recorded  C-V  characteristics  of  the  MOS  capacitor,  Fig.  1, 
and  from  the  normalized  minimum  capacitance  versus  oxide  thickness  for 
ideal  MIS  diodes  under  high  frequency  condition.  Fig.  2;  determination 
of  the  wafer  doping  after  gallium  diffusion  was  possible  [Ilf. 


Cmin 

^max 


118 

125 


0.944 


(1) 


The  oxide  thickness  can  be  determined  from 

A (2) 

Cmax  “ Cox  = eer  8ox 

where,  Cox  stands  for  oxide  capacitance;  A stands  for  area;  dox  stands 
for  oxide  thickness;  and  eer  stands  for  the  electric  permittivity. 

The  capacitor  had  a diameter  of  40  mils  and  an  area  of  (8.107  x 10  cm"), 

thus; 

2.83  x 10-15 

dox  = Cox  = 2260A°  (3) 

and  from  Fig.  2,  the  doping  is 


Na  * 10l8cm“3* 


(4) 


16  -3 

Since  the  original  substrate  doping  was  loss  than  10  cm  ; then 

jo  -3 

we  conclude  that  the  resultant  gallium  doping  is  around  10iOcm 


(Volts) 


d(A) 


Normalized  minimum  capacitance  versus 
oxide  thickness  for  ideal  MIS  diodes 
under  high-frequency  condition. 

AFTER  SZE  (12) 


FIGURE  2 


B.  GALLIUM  DIFFUSION  PROCESS  FOR  p-TYPE  WAFERS 
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The  following  is  a description  of  the  fabrication  procedure  of 
MOS  capacitors  and  MOSFET's  on  p-type  silicon  substrates. 

B.l  Fabrication  of  MOS  Capacitors  (Boron,  Gallium) 

The  following  is  a description  of  the  fabrication  procedure 
for  MOS  capacitors. 

1.  Oxide  growth  (4000A°  Si02>  - See  Fig.  (3). 

The  p-type  Silicon  Wafer  is  processed  as  follows  in  the 
oxidation  furnace. 

a)  20  Minutes  dry  oxidation  at  1100°c. 

b)  60  Minutes  wet  oxidation  at  1100°c. 

c)  20  Minutes  ary  oxidation  at  1100°c. 

d)  30  Minutes  aneal  in  Nj  at  1100°c. 

2.  Coat  with  Gallium  Silicafilra  - See  Fig.  (4). 

The  wafer  is  placed  on  the  spinner.  Eight  drops  of 
Gallium  Sillcafllm  are  depositied  on  the  wafer  while 
spinning  at  3000  rpm  for  ten  seconds. 

The  film  is  left  to  dry  at  room  temperature  for  one 
hour. 

3.  Gallium  Diffusion  - See  Fig.  (5)  and  (6). 

The  wafers  are  stacked  and  placed  in  the  diffusion 
furnace;  in  a pure  N2  ambient. 

a)  Predeposition — 30  minutes  at  1200°c. 

b)  Remove  the  Gallium  Sillcafllm  by  etching  in  a 
buffer  solution  (1:4  HF:NH4F)  for  one  to  one 
and  a half  minutes. 

c)  Drive  in  diffusion — 60  minutes  at  1150°c. 

4.  Metallzation  and  photolithography  to  obtain  a capacitor 
with  a guard  ring  for  surface  accumulation-See  Fig.  (7). 
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FIGURE  3 
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FIGURE  4 
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Third  Stc 


1200°  C 


Gal.  film 
Si02  — 


Si  Wafer 


Predeposition  at  1200°  C for 
30  minutes 


g ^ buffer  wafer 
L"  ' of  the  same  type 
JJ  as  the  coated  waffers 

Stacked  wafer  arrangement 
in  the  diffusion  furnace 


FIGURE  5 


Fourth  Stet 


II50°C 


Si  Wafer 


Drive  in  diffusion  for  60  min. 


FIGURE  6 
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FIGURE  7 
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B.2  Fabrication  of  Double  Doped,  MOS-FET  (Boron,  Gallium): 

The  following  is  a description  of  the  manufacturing  procedure 
for  MOSFET's. 

1.  Clean  surface  of  Si-p-type  wafer.  Boron  doped  - See  Fig.  (8). 

2.  Oxide  growth  7000A°  SIO2  - See  Fig. (9). 

a)  20  Minutes  dry  oxidation  1100°c. 

b)  180  Minutes  wet  oxidation  1100°c. 

c)  20  Minutes  dry  oxidation  1100°c. 

d)  30  Minutes  aneal  in  dry  No. 

3.  Photolithography: 

Photoresist  masking  step  to  open  the  source  and  drain 
diffusion  areas  in  the  oxide  - See  Fig.  (10). 

4.  Phosphorus  diffusion:  See  Fig.  (11). 

a)  Coat  the  surface  of  the  wafer  with  phosphoru.s  silica!  ilm. 
Deposit  5 drops  of  phosphorus  sllicafilm  on  the  surface 
of  the  wafer  while  spinning  at  3000  rpm  for  ten  seconds. 
(This  process  is  repeated  two  times).  The  wafer  is 

then  baked  at  200°c  for  15  minutes  to  dry  the  phosphorus 
film. 

b)  Predeposition:  one  hour  at  1100°c. 

c)  Remove  the  phosphorus  sllicafilm  by  HF  etch. 

d)  Drive  in  diffusion:  one  hour  at  1100°c  - See  Fig.  (12). 

e)  Etch  4000°A  of  the  oxide  without  exposing  the  gate 
area.  This  step  is  designed  to  protect  the  surface 
and  to  prevent  surface  Inversion  during  the  gate 
oxidation  process  due  to  phosphorus  evaporation. 

5.  Grow  3000° A gate  oxide:  Apply  the  same  procedure  outlined 
in  Step  2.  See  Fig.  (13). 

6.  Photolithography: 

Photoresist  masking  step  to  open  the  contact  holes  in 
the  oxide  of  the  source  and  drain  regions  - Sec  Fig.  (14). 
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FIGURE  15 


FIGURE  17 


7.  Gallium  Diffusion:  (see  Fig.  (15). 

a)  Deposit  8 drops  of  the  Gallium  silicafilm  on  the 
surface  of  the  wafer  while  spinning  at  3000  rpm 
for  ten  seconds. 

b)  Dry  the  gallium  film  at  room  temperature  for  one  hour. 

c)  Predeposition:  30  Minutes  at  1200°c. 

d)  Etch  the  silicafilm. 

e)  Drive  in  diffusion:  60  Minutes  at  1150°c  in  dry  Nj 
ambient . 

8.  Etch  700A°  of  S102  from  the  surface  to  remove  any  oxide 
formed  in  the  contact  holes. 

9.  Metalization  and  photolithography  to  establish  separate 
contacts  with  source,  drain  regions  and  gate  oxide. 

An  Improved  Process  for  the  Fabrication  of  Double  Doped  MOSFET 
The  MOSFET  fabrication  process  as  described  in  Section  B.2 
yields  devices  with  varying  gate  oxide  thickness  on  the 
same  wafer.  This  is  due  to  the  exposure  of  the  entire  wafer 
to  two  etching  steps  after  the  gate  oxide  growth.  However, 
by  combining  the  gate  oxide  growth  and  gallium  drive-in  within 
a step,  the  gate  oxide  uniformity  can  be  preserved.  In  this 
case  the  gate  oxide  growth  was  carried  out  as  follows: 

15  min.  dry  oxidation 
60  min.  wet  oxidation 
20  min.  anoal 

Sulfuric  acid  was  used  to  remove  any  remaining  gallium  at 
the  surface  of  the  wafer  after  Che  etching  step  of  the 
gallium  silicafilm.  This  was  found  to  enhance  considerably 
the  yield  of  the  process,  by  reducing  the  number  of  gate  shorts. 
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Therefore,  by  advancing  the  gallium  diffusion  process  to 
immediately  follow  the  phosphorus  diffusion,  by  using  sulfuric 
acid  wash  after  Ca  silicafilm  etch  and  by  combining  the 
gate  oxide  growth  and  gallium  drive-in  within  a step  a 
considerable  improvement  in  the  process  yield  and  in  the 
uniformity  of  device  characteristics  can  be  obtained. 

C.  Gallium  Doped  MOSFET  At  Room  Temperature 

The  MOSFET' s that  were  fabricated  according  to  the  process 

described  above,  were  mounted  and  measured  at  room  temperature. 

The  I-V  characteristics  of  one  of  them  is  shown  in  Fig.  (18). 

By  measuring  the  threshold  voltage  at  various  back  gate  voltages, 

17  -3 

the  Ga  concentration  was  estimated  to  be  = 10  cm 

Two  MOSFET  devices  were  selected  for  low  temperature  measurements, 
Device  A and  B.  The  dependence  of  threshold  voltage  on  back  gate 
bias  at  room  temperature  (Fig.  19)  enabled  us  to  deduce  the  total 

substrate  doping. 

The  slope  of  the  lines  in  Fig.  19  * 


from  equation  (5)  and  Fig.  19  we  determined  that  the  total  sub- 
strate doping  is  2.25  x 10^/cm^  for  device  A,  and  1.4  x lO^/cm-^ 
for  device  B. 
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III.  DEVICE  OPERATION 

Gallium  doping  introduces  an  acceptor  level  in  the  silicon  band  gap. 

The  energy  level  of  the  gallium  center  is  shown  in  Figure  20.  The  ioni- 
zation energy  E^  of  the  gallium  center  is  a measure  of  the  energy  required 
to  release  a hole  to  the  valence  energy  band. 

Gallium  is  an  impurity  from  group  III  of  the  periodic  table,  whicli 
requires  an  extra  electron  to  complete  the  valence  bonding  requirements 
[13].  The  gallium  center  receives  the  extra  electron  from  the  neighboring 
silicon  atoms  and  thereby  create  holes.  The  ionization  energy  of  the 
gallium  center  in  silicon  has  been  determined  by  Hall  measurements  [14,15, 
16,171.  In  this  case,  the  density  of  the  thermally  liberated  hole  which 
is  inversly  proportional  to  the  Hall  constant,  has  a temperature  dependance 

of  the  form  expC-E^/KT),  where  K is  the  Boltzmann's  constant  and  T is 

j 

the  absolute  temperature  [18].  The  gallium  ionization  energy  was  found 
to  be  0.065eV  [16].  The  gallium  ionization  energy  was  also  determined 
[19,20]  by  an  optical  method.  The  minimum  infrared  photon  energy  required 
to  cause  the  ionization  of  the  gallium  center  at  low  temperatures  was 
measured  [19]  and  found  to  be  0.071eV.  This  small  ionization  energy  for 
gallium  in  silicon  necessitates  very  low  temperatures  of  operation. 

Throughout  the  quoted  work  [14,15,16,17,18,19,20]  the  author  has 
found  no  reference  to  the  effect  of  the  electric  field  on  these  measured 
values  of  the  ionization  energies  for  gallium  in  silicon  or  for  any  ol 
the  other  shallow  acceptors  that  were  under  investigation.  Therefore, 
the  above  quoted  values  for  the  gallium  ionization  energy  should  not  be  com- 
pared to  the  values  we  calculated  fro  the  activation  energy  of  the  gallium 


centers . 


The  author  will  attempt  to  relate  the  above  quoted  E^  values  to 
values  we  measured  in  his  discussion  of  the  Poole  Frenkel  effect  in 
Section  V-A. 

Tne  gallium  center  in  silicon  at  low  temperatures,  has  two  distinctly 
different  charge  states.  The  most  positive  charge  state  corresponds  to 
the  neutral  center,  i.e.  the  center  when  occupied  by  a hole.  The  most 
negative  charge  state  corresponds  to  the  ionized  center  i.e.  the  center 
after  emitting  a hole  to  the  silicon  valence  band. 

The  ionization  of  the  gallium  center  can  be  either  due  to  thermal 
excitation,  or  to  incident  photons  with  energies  that  correspond  to  the 
far  infrared  radiation  wavelength  range.  The  effect  of  the  two  charge 
states  of  the  gallium  centers  on  the  operation  of  the  double  doped  gallium 
boron  MOSFET  can  only  be  observed  at  temperatures  below  50°K.  The  basic 
physical  phenomena  determining  the  operation  of  the  IRFET  can  be  described 
by  the  equations  of  the  Shockley-Read-Hall  model.  It  is  used  for  analyzing 
the  effect  of  the  two  charge  states  of  the  gallium  centers  in  the  depletion 
region  on  the  threshold  voltage  and  channel  conductance  of  the  MOSFET. 

The  gallium  doped  N-channel  MOSFET  (IRFET)  employs  gallium 
photoionization  to  modulate  the  conductance  of  the  channel  between  the 
source  and  drain.  A schematic  diagram  of  a cross  section  in  the  gallium 
doped  MOSFET  is  shown  on  Figure  21.  The  temperature  of  operation  will  be 
so  selected  such  that  the  gallium  centers  are  not  ionized. 


When  the  IRFET  is  used  in  an  imaging  array,  it  will  be  turncd-on  during 
the  period  prior  to  the  read  out.  This  period  will  be  referred  to  in 
our  discussion  as  the  read-in  period.  If  infrared  radiation  in  the  8 
to  14  micrometer  wavelength  range,  illuminates  the  array  during  the  read- 
in  period,  the  neutral  gallium  centers  in  the  space  charge  region  will  be 
photoionized  by  the  emission  of  a hole  to  the  valence  band.  Consequently, 
the  source  to  drain  resistance  of  the  MOSFET,  which  is  a function  of  the 
number  of  ionized  gallium  centers  in  the  space  charge  region,  will  be 
modulated  by  the  far  infrared  radiation.  The  source  to  drain  resistance 
of  the  IRFET  is  measured  during  the  read-out  period,  and  converted  to  a 
digital  signal  for  processing  and  display. 

The  gallium  doped  MOSFET  operates  consequently  as  an  integrating 
detector  during  the  read-in  period.  The  gallium  centers  will  be  partially 
or  totally  ionized  at  the  end  of  the  read- in  period  and  they  will  remain 
in  the  most  negative  charge  state  so  long  that  the  device  is  turned  on. 

The  read  in/out  cycle  can  then  be  repeated. 

The  physics  of  the  interaction  between  the  center  and  the  valence  band 
can  be  interpreted  in  terms  of  a hole  capture  process,  Cp,  an  optical 
emission  process,  e^,  and  a thermal  emission  process,  ep. 

The  definition  of  a depletion  region,  implies  the  nonexistence  of 
free  carriers  (holes  and  electrons).  Hence,  when  the  device  is  on,  the 
mechanism  of  hole  capture  for  the  gallium  centers  within  the  depletion 
region  can  be  ignored.  The  capture,  Cp  is  important  only  during  the 
reset  operation  of  the  device. 

At  sufficiently  low  temperatures,  the  optical  hole  emission  process 
from  the  gallium  center  to  the  valence  band,  e"  ,s  dominant.  The  thermal 
hole  emission  rate  will  be  effectively  eliminated  by  the  proper 
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selection  of  the  temperature  of  operation.  The  time  constant  (t),  of 
the  gallium  ionization  in  the  IRFET  will,  therefore,  be  equal  to  the 
reciprocal  of  the  optical  emission  rate,  e°. 

1 

T ■ e“  (6) 

A strong  field  dependence  has  been  observed  in  the  measurements  of 
the  thermal  emission  rate  for  holes  from  the  gallium  centers.  A detailed 
discussion  of  this  observation  can  be  found  in  Section  IV. 

The  temperature  of  operation  must  be  selected  low  enough  to  allow 
for  gallium  deionization.  But  it  must  not  again  be  too  low,  to  avoid 
the  deionization  of  the  shallow  level  acceptor  (boron)  in  the  substrate. 

The  ionization  energy  of  boron  lias  been  measured  [)6l  and  found  equal  to  0.045  eV 
The  complete  deionization  of  boron  will  raise  the  substrate  resistance 
and  may  disturb  the  accuracy  of  the  signal  processing. 

At  the  temperature  of  operation  the  thermal  time  constant  should 
be  at  least  two  orders  of  magnitude  longer  than  the  read-in/out  cycle. 
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IV.  KXPEK  I MKNTAL  SET-UP 
A.  Technique  ol'  Measurements 

The  MOSFET  devices  were  mounted  in  the  sample  holder  of  a cryogenic 
refrigerator,  which  is  basically  a liquid  hydrogen  cooled  cold-finger  in 
an  evacuated  metal  chamber. 

A diode  sensor  mounted  on  the  same  cold-finger  was  used  in  monitoring 
the  temperature  of  the  device.  A radiation  shield  cooled  to  80°K  had  an 
aperture  facing  the  filter.  The  filter  was  mounted  at  the  same  temperature 
to  limit  background  photon  flux. 

The  experimental  set  up  around  the  device  is  shown  in  Figure  22. 

The  device  was  triggered  by  a square  wave  on  the  gate.  A voltage  pro- 
portional to  the  source  to  drain  current  was  displayed  on  an 
oscilloscope.  The  waveforms  were  photographed  at  various  thermal, 
optical  and  electrical  bias  conditions.  The  photocurrent,  time  con- 
stant and  other  two-charge  state  characteristics  were  determined  from 
those  photographs.  This  technique  of  measurements  was  necessary  due 
to  the  fast  time  constant  of  gallium  ionization  at  the  temperatures  of  operation. 

All  voltages  were  monitored  on  FLUKE  DVM's.  The  MOSFET  and  the 
diode  sensor  were  biased  by  Hewlett  Packard  6216A  power  supply  units. 

The  gate  drive  voltage  was  obtained  from  the  bipolar  operational  power 
amplifier  (KEPCO) . The  waveforms  of  the  source  to  drain  current  were 
displayed  on  Tektronix  oscilloscope  typo  Shill  with  a differential 
amplifier  type  2A63. 

The  MOSFET's  used  in  this  experiment  had  circular  structures,  with 
outer,  middle  and  inner  rings  representing  source,  gate  and  drain 
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respectively.  The  dimension  of  those  rings  were  as  follows: 

The  source  outer  diameter  - 96  miL. 

The  source  inner  diameter  ~ 61  miL. 

The  gate  outer  diameter  a 66.5  miL. 

The  gate  inner  diameter  - 27.7  miL. 

The  diameter  of  the  drain  circle  » 28.5  miL.  Therefore  W-tt(  61+28 . 5) /2 

- 140  MiL,  L"(61-28.5)/2*16.25  and  W/L  = 8.17 
B.  Temperature  Measurements 

A diode  sensor  was  used  to  measure  the  temperature  of  the  sample 
holder  in  the  cryogenic  refrigerator.  This  diode  sensor  was 
calibrated  in  the  range  between  4 . 2°K  and  100°K . The  device  was 
purchased  from  LAKE  SHORE  CRYOTRON1CS,  INC.,  Model  No.  DT-200P/GR- 1 . 

The  diode  was  biased  in  the  forward  direction  by  10  microamperes. 

A digital  voltmeter  was  used  to  read  the  voltage  drop  across  the  diode. 

The  measured  voltage  drop  was  then  converted  through  the  manuf ac turer1  s 
calibration  table  into  temperature  readings. 

The  device  had  very  small  mass,  dimensions  and  dissipation  in 
the  range  of  20  microwatts, to  avoid  over  loading  the  system  thermally. 

We  have  neglected  in  our  temperature  measurements  the  error  resulting 
from  the  limited  input  impedance  of  the  DVM,  in  parallel  with  the  diode. 

We  have  also  assumed  that  the  temperature  at  the  corner  of  the  sample 
holder,  where  the  diode  was  mounted,  is  identical  to  the  temperature 
of  the  device.  This  assumption  is  practically  justified,  since  the 
sample  holder  was  designed  with  minimum  dimensions.  The  diode  sensor 
measured  a combination  of  both  the  lead  and  case  temperatures.  The 
lead  winding  around  the  frame  had  to  be  done  with  extreme  care  to  insure 
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accurate  temperature  readings.  The  manufacturer  's  i .i  I ibr.u  ion  t . ■ | w.is 
checked  at  77°K  and  300°K  before  the  device  was  used. 

Before  the  diode  sensor  was  purchased,  we  attempted  to  measure 
the  temperature  by  a commercial  silicon  diode.  The  same  10  micro- 
amperes forward  bias  arrangement  was  used.  Calibration  of  the 
voltage  drop  versus  temperature  was  made  at  300°K  and  77°K.  These 
two  calibration  points  were  used  to  develop  a temperature  correction 
factor.  The  temperatures  were  calculated  from  the  ideal  1-V  diode 
equation  under  forward  bias  condition.  This  temperature  measurement 
technique, using  a commercial  diode, left  a lot  to  be  desired  in  terms 
of  accuracy,  which  ultimately  caused  us  to  decide  in  favor  of  using 
the  LAKE  SHORE  diode  sensor. 

C . Circular  Variable  Filter  [21] 

Multilayer  coating  design  allowed  the  production  of  any  filter 
in  the  "CVF"  configuration.  Optical  interference  coatings  are  vacuum 
deposited  on  circular  germanium  substrates.  Film  thickness,  and 
therefore  wavelength,  varies  linearly  with  angular  position. 

The  Infrared  circular  variable  filter  was  purchased  from  OCl I 
"Optical  Coating  Laboratory  Inc."  The  infrared  range  is  covered  by 
a set  of  three  90°  segments.  The  wavelength  coverage  cxLends  from  2.5 
micrometer  up  to  14.5  micrometer.  In  the  experimental  work  done  on 
the  gallium-boron  doped  1RFET;  only  segment  three  was  used.  Tills 
segment  had  the  following  specifications:  Wavelength  range,  X«7.9Pm  to 
14.5  pm,  half  bandwidth  < 1.82,  transmittance  in  band  > 302  peak,  and 
blocking  for  wavelengths  below  the  bandpass;  corresponds  to  a 
transmission  £ 0.12.  This  value  is  also  applicable  for  wavelengths  up  to 
17.0  pm,  when  a blocking  element  is  used. 
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Employing  CVF's  eliminates  the  need  for  multiple  filters,  and 
allows  for  rapid  scan  capabilities  with  high  resolution.  More  important, 
those  filters  are  capable  of  being  operated  at  cryogenic  temperatures, 
down  to  4°K.  The  problem  of  1RFET  discharge  by  photon  flux  I roiti 
filters  located  in  300°K  ambient,  can  be  eliminated  if  the  filter  is 
attached  to  the  77°K  radiation  shield.  A larger  outer  chamber  for  the 
cryogenic  refrigerator  had  to  be  designed  to  accomodate  the  filter, 
and  allow  for  segment  rotation. 

The  manufacturer  quotes  wavelength  variations  with  temperatures 
of  an  order  less  than  1.0%  of  the  center  wavelength  of  the  bandpass, per 
100°C.  Center  wavelength  has  a tendency  to  shift  to  shorter  wavelengths 
with  decrease  in  temperatures  [21  ].  The  magnitude  of  the  wavelength 

shift  is  dependent  on  the  variation  in  the  refractive  indices  of  the 
material  used  as  a result  of  temperature  variation.  The  manufacturer 
also  guarantees  the  linearity  of  the  transmittance  band  vs.  angle, within 
3%  range  for  points  at  the  same  radius. 

D.  Cryogenic  Refrigerator 

Model  AC-2  CRY0-T1P  refrigerator  was  used  to  carry  out  the  low 
temperature  measurements  on  the  Ca-B  doped  MOSFET.  The  system  is  made 
by  Air  products  and  Chemicals  Inc.  Model  AC-2  is  a minuature  open  cycle 
Joule -Thomson  refrigerator  with  coil  finned  heat  exchanger.  The  net 
refrigeration  capacity  Is  4 watts  maximum  at  23°K  122],  the  specified  tem- 
perature range  Is  16°K  to  300°K,  and  the  refrigeration  is  achieved  from 
standard  compressed  cylinder  gasses  (Nitrogen  and  Hydrogen). 

The  cooldown  time  was  45  minutes  on  the  average.  The  minimum  tem- 
perature achieved  with  the  WMX-1A  vacuum  shroud  and  radiation  shield  was 
21°K.  The  radiation  shield  is  kept  at  approximately  the  liquid  nitrogen 
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temperature  (23),  thus  providing  a shield  at  about  80°K  to  the  liquid 
hydrogen  (20.4°K)  cooled  sample  holder.  The  maximum  refrigeration  capacity 
is  achieved  with  this  configuration,  and  any  increase  in  the  size  or  mass 
of  the  vacuum  shroud,  radiation  shield  or  specimen  holder  is  bound  to  in- 
crease the  lowest  achievable  temperature  [24]. 

The  optical  measurements  on  the  Ga-B  doped  MOSFET  had  to  be  con- 
ducted with  a cold  filter.  A larger  vacuum  shroud  was  designed  and  pro- 
duced in  the  U of  A workshop,  to  enable  us  to  mount  and  move  the  circular 
variable  filter  in  the  80°K  temperature  station.  The  operation  experience 
with  the  AC-2  Cryo-Tip  and  the  heat  exchanger  model  No.  AC-2-110,  showed 
that  the  running  time  per  220  Cu.  ft.  N2  cylinder  is  2 hours  including 
the  initial  load  of  cooldown.  The  running  time  per  220  cu.  ft.  cylinder 
is  3 hours.  The  and  1^  cylinders  were  charged  to  2400  PSI  and  2100  PSl 
respectively. 

The  specimen  holders  were  designed  with  minimum  dimensions,  to  re- 
duce the  thermal  load  on  the  cold  finger  and  to  obtain  the  minimum  possible 
coldf inger  temperature.  The  specimen  holders  were  fabricated  from  copper 
and  oxygen  free  copper  materials  due  to  their  good  thermal  conductivity 
at  low  temperatures.  Vacuum  of  4X10  ^ torr  is  necessary  to  limit  heat 
leak,  and  obtain  the  lowest  operation  temperatures. 
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V.  RESULTS  AND  ANALYSIS 

A.  Thermal  Emission  Rate  and  Poole-Krenklc  Effect 

To  fully  characterize  the  device,  one  should  determine  the  range 
of  temperatures  where  the  ionization  of  the  gallium  center  does  not 
occur  instantaneously. 

The  device  was  triggered  by  a square  wave  on  the  gate  with 
Vgg  - + 40V,  Vgc  « 0 and  VDS  * +10  volts.  A voltage  proportional  to  the 
source  to  drain  current  was  displayed  on  the  oscilloscope.  The  decay 
of  source  to  drain  current,  with  a time  constant,  from  the  initial  value 
to  its  final  value  was  photograohed  at  various  temperatures.  The  device 
was  fully  screened  from  any  radiation  source  during  this  process. 

The  magnitude  of  the  decay  in  the  source  to  drain  current, 
corresponds  to  the  limit  value  of  the  photocurrent , Al^g  [3].  Wo  have 
associated  the  current  decay  rate  with  an  exponential  function 
and  the  time  constant  was  determined  for  each  temperature  and  plotted 
in  Figure 23  versus  1000/T. 

Since  the  thermal  emission  process  is  the  only  ionization  process 
producing  this  decay,  we  concluded  that  the  time  constant  is  equal  to 
the  thermal  time  constant  and  the  reciprocal  of  the  thermal  emission 

rate. 

1 
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In  reality,  the  decay  corresponded  more  to  a combination  of  a fast  ex- 
ponential function  superimposed  on  a slower  exponential  function,  partic- 
ularly at  very  low  temperatures.  The  time  constants  of  both  of  them 
were  temperature  dependent,  with  the  fast  time  constant  virtually  dis- 
appearing at  higher  temperatures. 


rjrrarrg 


TIME  CONSTANT  FOR  THERMAL  EMISSION 
r • |/«L  (SEC) 
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T { K°) 


1000/  T 
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The  fast  time  constant  was  interpreted  as  being  associated  with 
centers  which  are  physically  near  to  the  drain  region.  Where  .as  the 
slower  time  constant  is  for  the  gallium  centers  in  the  rest  of  the  de- 
pletion regio  ..  However,  this  interpretation  does  not  exclude  the 
possibility  that  some  of  the  fast  decay  might  have  been  partly  due  to 
boron  deionization  or  other  causes. 

The  thermal  time  constant  can  be  described  by  the  following  equation. 

_ 1 . +AE/KT 

T - — — * A e 


where  AE  is  the  activation  energy 
constant  and  T is  the  temperature 
can  be  determined. 

For  the  lightly  doped  device 


at  1000  = 33 
T 

at  1000  » 46 
T 


therefore  ■ 1 * A exp 

-2 

10  » A exp 


for  the  center,  K is  the  Boltzmann 
in  °K.  From  Figure  23  , A and  AE 


(A):  N - 2.25  x 10lf>cm~3 
-2 

T ^ “ 10  sec. 


T * 1 second 


AE  46 
K 1000 


AE  33 
K 1000 


(9) 

(10) 


hence  by  simple  division,  we  get 


or 


therefore 

and 


100  - 


exp 


AE  1 3 

k i ooo 


4.61 


AE  13 
K 1000 


AE  - 0.0307  eV 
A - exp  (-354.24  (H-  )) 

nooo" 

* 8.4  x 10  seconds 


(ID 


(12) 


(13) 

(14) 
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which  is  the  limit  value  of  the  time  constant  at  very  high  temperatures, 

17  -3 

Similarly  for  the  heavily  doped  device  ( B) : * 1.4  x 10  cm 


therefore 


at 

1000 

- 29 

r = 

f- 

10  sec 

T 

L 

at 

1000 

* 43.4 

Tt  * 

10  3 sec . 

T 

T 

t 

•-H 

1 

o 

H 

N 

= A exp 

(— 
v K 

43.4 
1000  } 

(15) 

io-3 

« A exp 

.AE 
{ K 

29  ) 
1000' 

(16) 

and  by  simple  division,  we  get 

inn  ,AE  14.4, 

100  - exp  ( K 100Q) 


Therefore 


, ,,  AE  14.4 
4‘61  ‘ ~K  1000 


AE  - 0.02772  eV 

a 

t n nn  i a ° 


and  A * 9.38  x 10  seconds  (20) 

which  is  the  limit  value  of  the  time  constant  at  very  high  temperatures. 

The  thermal  time  constant  was  found  to  be  strongly  field  dependent. 

The  decay  was  several  times  faster  with  a very  small  back  gate  bias  po- 
tential applied. 

The  ionization  energy  of  gallium  was  measured  by  others  and  found 
equal  to  0.071  eV. 

In  the  cases  of  double  doped  MOSFET  employing  gold  or  indium  as 
second  dopant  [ 2],  [ 5],  the  activation  energy  determined  by 
the  thermal  time  constant  and  the  temperature;  was  found  to  be  almost 
identical  with  the  ionization  energies  of  gold  and  indium  0.35  eV  and  0.155  eV 


f ; 
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respectively.  However,  the  deviation  observed  here,  in  the  case  of 
the  gallium  center  is  explainable  by  the  Poole -Frenkel  effect.  In 
other  words  0.072  eV  - 0.0307  eV  - 0.0413  eV  or  57. 4%  out  of  the 
total  ionization  energy  of  the  center  was  supplied  by  the  electric 
field  across  the  depletion  region,  rather  than  by  the  thermal  energy. 
Simply,  since  the  gallium  center  is  a shallower  acceptor,  it  is  more 
sensitive  to  the  electric  field  contribution  during  the  ionization.  The 
measurements  were  repeated  with  different  source  to  drain  voltage  V 

I/O 

and  consistent  results  were  obtained. 

The  rate  of  hole  emission  will  not  only  be  dependent  on  the  location 
of  the  impurity  energy  level  w.r.t.  the  valence  band,  but  also  on  the 
density  of  unoccuppied  states  within  the  valence  band  [10].  Therefore, 
we  have  attempted  to  match  the  thermal  responce  data  of  Fig.  23,  to  the 
following  equation: 

Tt  - ~ - A T_3/2  exp  (+AE)/kT.  (21) 

e 

P 


This  equation  should  be  valid  under  the  assumption  of  constant  thermal 
velocity  in  the  temperature  range  under  investigation  by  ignoring  the 
carrier  "heating"  resulting  from  high  electric  field. 

(vth  - /3KT / m - 3.65X106cm/sec  (at  29.4°K)  (22) 


at  29.4°K) . 

For  the  lightly  doped  device  (A):  Nj  * 2.25X10^cm  ^ we  obtain: 


and 


&E  - 0.02725eV 
A - 4.94X10-5  sec.  °K-3/2 


(23) 

(24) 
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Similarly  for  the  heavily  doped  device 


(B):  N - 1.4X1017cm-3 


we  obtain: 


AE  *■  0.023946eV 


A - 6.47X10_3sec.°K_3/2 


(25) 


In  the  following  section  the  field  depend. ince  of  the  emission 
coefficient  for  carriers  from  shallow  acceptor  levels  in  the  depletion 
region  of  Si-MOSFET  behind  a strongly  inverted  surface  will  be  derived. 

The  time  dependance  of  the  gallium  charge  state  in  the  depletion 
region  of  the  IRFET,  can  be  determined  from  the  Shockley-Read-Ha 1 1 (SRH) 
model.  The  total  doping  concentration  will  be  defined  as 


NTT  - NA  + Nca 


(26) 


in  the  steady  state  inversion  condition.  Where  NA  is  the  concentration 
of  the  boron  atoms,  and  Nq3  is  the  concentration  of  the  gallium  atoms. 
In  the  transient  condition  that  follows  the  application  of  an  inversion 
voltage  to  the  gate  plate. 


N™  - N_  + 

TT  T T 


(27) 


where  N.p  is  the  concentration  of  the  negative  or  ionized  acceptor  centers, 
and  PT  is  the  concentration  of  the  centers  in  the  neutral  charge  state. 

W.  Shockley  and  W.  T.  Read,  lr.  have  studied  the  statistics  of  the 
recombination  occurring  through  the  mechanism  of  trapping  [25],  for  both 
equilibrium  and  nonequilibrium  steady  state  conditions.  They  also  analyzed 
the  rate  of  generation  in  the  space  charge  regions.  C.  T.  Sah  [26]  derived 
the  equivalent  circuit  model  for  the  SRH  center.  The  equivalent  circuit 


I w. 


y»'  ' 


mentioned  above  is  a powerful  tool  in  determining  the  response  under 
transient  conditions  superimposed  on  an  arbitrary  steady  state  force. 

The  rate  equations  of  the  above  model  are  now  applied  to  the  transient 
response  of  the  shallow  acceptor-gallium  center  in  the  depletion  region 
of  a double  doped  gallium-boron  MOSFET,  which  operates  at  a specific 
low  temperature  region  which  would  allow  the  gallium  to  be  neutral  while 
maintaining  total  boron  ionization.  The  rate  equation  describing  the  change 
in  N with  respect  to  time  125,  26]  is  given  by 


T = C nP  - C pN  + e P - e N (28) 

nr  P 1 P T n 1 

where  e is  the  emission  rate  of  electrons,  e is  the  emission  rate  of  holes, 
n P 

is  the  capture  rate  of  electrons,  is  the  capture  rate  of  holes,  n is 
the  number  of  available  electrons  and  p is  the  number  of  available  holes. 

In  a depletion  region,  though,  there  are  no  free  carriers  available  for 
capture,  i.e.  n » p = 0.  Thus  the  rate  equation  becomes 


“dt  -epPT-enNT  <29> 

Our  initial  conditions  correspond  to  the  end  of  the  accumulation  period 
at  time  t*0.  During  this  accumulation  period  each  gallium  center  captures 
a hole.  Thus,  all  the  gallium  centers  enter  the  transient  period  in  their 
neutral  charge  state  (most  positive  charge  state),  where  as  all  the  boron 
centers  are  ionized  (negative  charge  state).  Therefore  at  t-0 

Nt(0)  - Na  (30) 

PT(°)  - NCa  (11) 

The  final  steady  state  condition  with  an  inverted  surface  should  be 


. »s  ♦ «Ca 

rT(fli.al)  > 0 


a***,-*:  ■ewnattWe'***  -iiSSiSttv 
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From  equation  (27)  and  (29);  we  obtain 


dN_ 


- e (N„_  - N._)  - e N_ 
dt  pv  TT  T n T 


(36) 


therefore 

dN_ 


dt  " 6pNTT  " (en  + fep)  NT 


(35) 


from  equation  (35),  we  obtain 
dN 


dt 


e N ~ Me+*) 

p TT  T p n 


(3b) 


and 


dN, 


N - p M 
T — N__ 
e +e  TT 
p n 


-(e  + e )dt 

n p 


(37) 


therefore 


Nt  - ep  _ Ntt  - C exp  (-  (ep  + en)t  ) 


(38) 


e +e 
p n 


where  C is  a constant. 

Equation  (38)  can  also  be  written  in  the  following  form. 


N_(t)  - C exp  - t/x  + 


N 


where 


e +e 
P n 


TT 


(39) 


e + e 
P n 


(60) 


We  evaluate  the  constant  C from  the  initial  conditions:  at  t-0  N » N 

T A 


therefore 


N - C + N 

A e +e  TT 

P n 


(61) 


and  hence 


C - N.  — N__ 

A e +e  TT 
P n 


e +e 
P n 
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V1  - \ A. 

e +e  p n 

P " 


Therefore  equation  (39)  can  be  written  ns  follows 


NT(t)  - (Na  %_  - NCa 

e +e  e +e 

p n p n 


] exp  -t/x  + ep  (Na  + NCa> 


e +e 
P » 


In  the  final  steady  state  condition 

KT(fl™l>  - (»  + « ) <«> 

p n 

But  since  the  gallium  acceptor  level  is  much  closer  to  the  valence  band 

than  to  the  conduction  band,  then  e >>  c and  equations  (43)  and  (44) 

p n 

yield  the  following  two  equations: 

*!<*>  ’ - Nc,  r'/T  + \ * »c.a  <“> 


NT(final)  a NA  + NCa  (46) 

where  the  average  time  constant  T ■ — under  nonequilibrium  conditions. 

Pp  t r 

e is  comprised  of  both  the  thermal  field  assisted  emission  rate  c ’ and 
P P 

the  optical  emission  rate  e”,  therefore 

e - i - e t,e  + e°  (1/sec)  (47) 

PIP  P 

substituting  for  both  the  thermal  field  assisted  emission  rate  and 

the  optical  emission  rate  in  equation  (47)  yielcfethe  equation  of  the  general 
emission  rate  of  an  electric  field  tunnable  infrared  detector. 

« -A.-'V"  ♦Be+flE"/"T  «l) 

P P 

where  A,B  have  the  dimensions  of  — 

sec 

A + B : is  the  emission  rate  at  low  electric  field  and  high 

temperature  while  the  device  is  operating  in  the  dark 

o 2 

Op  : i 8 the  photoionization  cross  section  (cm  ) 

2 

4>  : photon  flux  0/cm  . sec 


• V «*• 


U2 


AE 

t 

E 

n 


B 


the  thermal  activation  energy 
the  average  or  mean  electric  field 

a constant.  The  value  of  this  constant  depends  on  the 
degree  of  nonlinearity  of  the  electric  field  in  the 
depletion  region  of  the  MOSFET. 
the  Poole-Frenkel  constant 

6 - (q3/irK8lCo)'S  (49) 


K : the  Boltzmann  constant 

T : the  temperature  in  °K 

As  shown  earlier;  we  measured  the  response  of  the  gallium  center  in 
the  depletion  region  of  a MOSFET  behind  a strongly  inverted  surface.  If  the 
thermal  activation  energy  AE^  determined  from  these  measurements  were 
added  to  the  barrier  lowering  due  to  electric  field,  the  result  would  be 
related  to  the  ionization  energy  of  the  gallium  center,  [16,  19]. 

The  Poole-Frenkel  theory  predicts  that  the  electric  field  will  lower 

u 

the  coulomb  barrier  by  an  amount  proportional  to  E , thus  increasing  the 
probability  of  thermally  induced-emission  [8] . By  applying  a simple  one 
dimensional  Poole-Frenkel  model  we  will  illustrate  how  the  probability 
of  thermally  induced  emission  is  enhanced  by  electric  field.  The  Poole- 
Frenkel  effect  is  very  similar  to  the  field  assisted  Schottky  emission 
f-8,9],  but  it  is  twice  as  strong  because  the  interaction  of  the  escaping 
carrier  is  with  a fixed  point  charge  rather  than  with  an  image  charge  |8|. 

Since  the  maximum  electric  field  across  the  depletion  region  in  both 
device  A and  B is  less  than  6 x 10"*  V/cm,  we  do  not  expert’  that  ionization 
by  carrier  tunneling  has  played  any  significant  role  in  the  device  operation. 
However,  the  thermal  velocity  can  be  much  larger  [8]  than  that  corresponding 
to  the  lattice  temperature  because  of  carrier  "heating"  by  electric  field 
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larger  than  100  V/cm  at  temperatures  below  77°K  |81. 

In  a n-channel  MOSFET  the  band  bending  at  the  surface  when  the 

device  Is  on  deve lopes  the  surface  potential,  4>  . The  surface  potential 

s 

is  approximately  equal  to 


*s  * 2tf  + VBG  (50) 

under  back  gate  bias  condition.  The  surface  potential  however,  is  also 
effected  by  the  drain  to  source  voltage.  If  the  later  results  in  a uniform 
spatial  electric  field  across  the  surface,  then  we  can  approximate  the 
surface  potential  as  follows: 


*8  - 2*f  + VBG  + 2 


DS 


(51) 


at  a point  midway  between  the  source  and  drain.  This  approximation  does 
not  take  in  account  the  effect  of  channel  pinch  off  on  the  linearity  of 
the  voltage  distribution  along  the  channel. 

The  maximum  depletion  layer  depth,  as  defined  in  the  MOS  theory  (27], 
would  be  equal  to 

Xn  - /2K  (£  O'"  (52) 

I'  qNT(0 

where  <f>g  is  the  surface  potential  as  defined  by  either  equation(50)  or  (51), 

N^Ct)  is  the  concentration  of  the  ionized  acceptors  as  defined  by  either 

equation  (43)  or  (45),  k * 12,  e * 8.85  x 10  ^ F/cm  and  o is  the 

si  o 

-19 

electronic  charge  * 1.6  x 10  coulomb.  We  conclude  from  equation  (52) 
that  the  depth  of  the  depletion  layer  is  a function  of  time  during  the 
transient  period  or  the  read-in  period.  The  maximum  depletion  layer  depth 
occurs  when  N^.(t)  **  N^,  and  the  minimum  occurs  at  the  end  of  the  gallium 
ionization  process  when  N^,(t)  ■ If  the  effect  of  space  charge 

on  the  voltage  distribution  across  the  depletion  layer  is  Ignored,  in 
other  words,  if  we  assume  that  the  voltage  drop  is  linearly  distributed 
in  the  direction  normal  to  the  silicon  surface,  then  the  electric  field 


ia  constant  along  the  depletion  layer  depth  and  equals 

* - !» 


(53) 


We  conclude  from  equation  (53)  that  the  electric  field  across  the 
depletion  layer  is  also  a function  of  time  as  a result  of  the  l line 
varying  depletion  layer  depth.  This  ultimately  means  that  the  emission 
rate  as  defined  in  equation  (48)  is  also  a function  of  time  through 
its  exponential  dependance  on  the  electric  field,  with  the  minimum 
emission  rate  occurring  at  the  beginning  of  the  ionization  process  and 
the  maximum  emission  rate  occuring  at  the  end  of  that  process.  This 
can  also  explain  why  did  the  decay  time  constant  appear  as  an  initial 
fast  decay  superimposed  on  a slower  decay. 

From  equation  (45)  we  will  calculate  now  the  mean  doping  in  the 

depletion  region  during  the  transient  under  the  assumption  that  e >>  e . 

<N  > - N (t) 

T J dt  - N A + 0.755  H,  (54) 


*0*  ‘ 


2|SiCo^ 


/ ^A-  0*755 


Now,  for  device  A,  at  low  temperature,  under  the  assumption  that  the 
fermi  level  lies  above  the  boron  energy  level  and  below  the  gallium  energy 
level  in  the  bulk  silicon  such  that  all  the  gallium  atoms  are  neutral 
except  when  the  device  is  turned  on. 

<Nt>  - 1.8  x 1016  + 0.755  (0.45  x 1016) 

■ 2.1  x 10^6  /cm"*  (57) 

and  using  the  surface  potential  approximation  of  equations  (50)  and  (51) 
we  calculate  the  mean  depletion  layer  depth  where  2<t>f  = 1 at  low  temperatures 
From  equation  (50)and  (55)  we  obtain 


vV"  09- 
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<v 


/ 2 X 12  x 8.85  k 10  l4(l  + v ) 


/- 


6 x 10  ^ x 2.1  x 10*° 


2.5  x 10-5  (1  + \'BG)’J 


(58) 


where  as  from  equation  (51)  and  (55)  we  obtain 


<XD>  - 2.5  x 10'5  (1  + VBG+  VDS)b 
from  the  equations  (56)  and  (58),  we  do  (luce  that 
<£>  > 3.98  x 104(1  + VBG) -1 


(59) 


(60) 


and  from  the  equations  (56)  and  (59)  , we  obtain 
<E>  - 3,98  x 104  (1  + VBG  + VDS)'5 


(61) 


Now  let  us  calculate  the  Poole-Frenkel  constant  B: 


e - 


uK  , £ 
si  o 


1.6  x 10 


-19 


3.14  x 1 2x  8.85  x 10 


i, 

1 2 U 1 

-14;  = 2.2  x 10  cV(cm)  2 


(62) 


KT  at  26°K 


1.38  x 10 


-23 


1.6  x 10 


-19 


26  = 2.24  x 10-ieV 


(63) 


KT  at  34. 7°K»  8.617  x 10  3(34.7°)  = 2.99  x 10~3t«V 


(64) 


from  (62)  and  (63) 


P 2.2  x IQ'*  = 9.82  x 10-2 

KT  J " 


2.24  x 10 


-3 


at  26°K 


from  (62)  and  (64) 


(65) 


P , 2.2  x 10 

] “ 


-4 


KT  J „ „„  . -3 

2.99  x 10 


7.36  x 10 


-2 


(66) 


at  34 . 7 K 
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Therefore  using  equation  (60)  and  both  equations  (65)  and  (66),  wo  coni  ludo 
that  the  rate  of  hole  emission  is  proportional  to  the  following 
e <»  exp  7T  E 4 

p M 

assuming  a constant  electric  field  across  the  depletion  layer. 

At  26°K  for  device  A we  get 

e « exp  9.82  x 10-2  x 199.5  (1  + V_„)!'‘  (6H) 

p r BO 

“ exp  19.6  (i  + V^)*4 
At  34.7°K  for  device  A wc  get 

e “ exp  14.6  (1  + V ) 4 (69) 

P DVj 

Equations  (68)  and  (69)  are  plotted  together  with  the  measured  data  on 
Figure  24.  The  slopes  of  the  resulting  lines  on  a semi-log  plot  match 
the  measured  points  very  closely. 

The  barrier  lowering,  at  (V  = 0),  can  be  calculated  from  the  Poole- 

BO 

F.  nkel  model 

A<t>pF  = 3^  = 2.2  x 10“4  x 199.5  = 0.0439eV 
Now  summing  the  barrier  lowering  with  the  measured  activation  energy 
we  get  the  ionization  energy  of  the  gallium  center. 

■ 0.02725  +0.02.39  - 0.071eV  (70) 

which  is  identical  to  the  ionization  energy  of  the  gallium  center  measured 
by  E.  flurstein  in  1955  (19]. 

Similarly  using  equation  (61)  and  both  equations  (65)  and  (66),  wo  obtain 
the  proportionality  between  the  hole  emission  rate  and  the  electric  field. 

At  26°K  for  device  A we  get 

VDS  k 

ep  « exp  19.6  (1  + V^,  + -!p)  * (71) 

and  at  34.7°K  for  device  A we  get 

V w 

exp  14.6  (1  + VBG  + )* 


e <*■ 
P 


(72) 


THERMAL  EMISSION  RATE,  en(sec) 


RG.24  POOLE  - FRENKLE  EFFECT 
v HELD  ENHANCED  EMISSION 


f 
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Equations  (71)  and  (72)  are  plotted  together  with  our  measurement. s on 
Figure  25.  The  slopes  of  the  resulting  lines  on  a semi-log  graph  make 
perfect  match  with  the  measured  points. 

The  barrier  lowering  at  V =0  and  V =8.9  v0^tj  Can  be  calculated 

Ij\r  Ui) 

from  the  Poole-*-Frenkcl  model 

i 8.9  , 

A<J>pf  - 0/e  - 2.2  X 10  X 199.5  x (1  + 2 V 

= 0.067  eV 


(73) 


The  ionization  energy  in  this  case  should  be 


AEt  + 


A<(> 


PF 


= 0.02725  + 0.067  = 0.094  eV 


(74) 


which  is  appreciably  different  than  the  previously  measured  ionization 
energies  for  the  gallium  center.  The  deviation  is  due  to  our  earlier 
assumption  that  is  linearly  dropped  across  a pinched  off  channel  in 
the  saturation  region. 

In  all  our  analysis  so  far,  we  have  ignored  the  space  charge  effect 
on  the  surface  potential.  In  other  words  we  assumed  that  the  electric 
field  is  constant  along  the  depth  of  the  depletior  layer.  In  reality 
however,  the  situation  is  more,  complex.  The  surface  potential  distribution 
is  obtained  from  the  solution  of  Poisson's  equation  in  the  space  charge 
region  [12]. 


K c 
s o 


3 

3x^ 


= -q  <Nt> 


which  have  a solution  of  the  following  form 


(75) 


♦ • ♦.  » - ^>i2  <«> 

which  yields  a potential  energy  distribution  of  the  form 

(P.E.)  - -q  ♦ [ 1 - l2  (77) 

» u 

and  an  electric  field  distribution  of  the  form 


(78) 


'><) 


On  the  other  hand  the  force  of  attraction  between  the  ion  and  the 

escaping  carrier  is  of  the  form 
2 

F - 

a o 


4"k  x 
si  e 


(79) 


The  work  done  by  the  carrier  In  the  course  of  its  transfer  from 
infinity  to  the  point  x,  corresponds  to  the  potential  energy  of  the 

carrier  at  point  x and  is  given  by 

2 

q 

(PE)c  - 

4tk  ,e  x 
si  o 


(80) 


(81) 


Therefore  the  total  potential  energy  as  a function  of  distance  is  given 
by  equation  (81).  Where  X=0  at  Si  - SiO^  interface. 

PE(x)  * 9 + q ^-s  * x) 

4mCsiEoI(<V  ‘ X)  Yl  <V2 

where: 

y <0Cjj  - x and  0 < x < <3tp>.  as  illustrated  on  Fgiure  (26-A).  Where  y Fig- 
ure 26  is  the  separation  between  the  ion  and  the  edge  of  the  depletion 
layer.  Therefore  the  location  of  the  lowering  x^  is  given  by  the  condition 


d[PE(x) ]/d  [<XD>  - x]  * o 

.» 


X_  - X 
D m 


X (at  y-0) 
m 


[[4TT<NT>j] 


1/3 


cm 


therefore  the  resulting  barrier  lowering  is  equal  to 

_q 

KsiCo 


| A4> (at  y“0)  | “ “ j 


> 1/3 

[—41 

16tt^  volts 


(82) 


(83) 


for  y i1  0 


<X  > - X - X-  X (0)  [ 1 - -£-]  3 

u mm  a 

m 


and 


A4>  ■ 


>!i  + "V  *.2  1 

Si  O IT 


(84) 


(85) 


In  other  words  both  the  position  and  value  of  the  lowering  are  functions 
of  the  ion  separation  from  the  depletion  layer  edge.  The  complexity 


CROSS  SECTION  OF  MOS- 
TRANSISTOR  OPERATING  IN  SATURATION 


of  the  real  situation,  makes  us  favor  the  use  of  the  simple  Poole- 
Frenkel  model,  but  barring  in  mind  that  there  is  some  error  evolved. 

D.  F.  Bentchkowsky  and  A.  S-  Grove  studied  electric  field  dis- 
turbance in  the  drain  to  channel  region  [28]  for  a device  operating  in 
saturation.  They  also  tried  to  avoid  solving  Poisson's  equation  is 
this  region  to  estimate  precisely  the  effect  of  the  two  transverse  field 
contributions  on  the  channel  conductance.  They  favored  the  use  of  a 
simple  less  accurate  model  over  the  complexity  of  the  solution  of  Poisson's 
equation  in  the  channel  region. 

The  solutions  of  the  equations  (84)  and  (85)  are  summarized  in  the 
following  table  and  plotted  on  figure  (26-B). 


y in  A° 

0 

20 

40 

60 

100 

2000 

2400 

xm  in  A° 

156 

170 

183.8 

197.8 

228.8 

2043.1 

2439.4 

A4>in  volts 

0.01156 

0.0126 

0.0137 

0.0148 

0.0175 

0.688 

0.972 

This  solution  was  obtained  by  substituting  < Ny(t)>  = 2.1  x lO^cm-^ 
for  device  A at  VgG  - 0.  The  most  interesting  tiling  about  those  cal- 
culations, is  that  they  suggest  that  some  of  the  gallium  concentration 
in  the  depletion  region  may  be  permanently  ionized  due  to  the  surface 
potential.  Naturally  those  are  the  gallium  atoms  nearest  to  the  surface. 


FIG.  26 B Graphical  Representation  of  the  Solution  of 
Roissans  Equation  in  the  Depletion  Region  of  Ga-B  doped  MOSFET 
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B.  Two  Charge  State  Device  Characteristics 

The  gallium-boron  MOSFET  device  A,  Nj  * 2.25  x lO^cm-^ 
was  tested  on  the  curve  tracer  (Type  576,  Tektronix  Inc.).  The  1-V 
characteristics  of  the  device  are  shown  in  Figure  27  . The  characteristics 
were  photographed  from  the  tracer  display  with  a time  exposure. 

The  dark  lines  correspond  to  the  initial  condition,  with  the  gallium  centers 
in  their  neutral  state,  immediately  after  the  surface  has  been  accumulated. 
The  dashed  lines  correspond  to  the  final  condition,  with  the  gallium 
centers  ionized  or  in  their  negative  charge  state.  This  decay  in  the 
I-V  characteristics  was  measured  at  24°K  and  without  any  back  gate  bias 
potential  applied.  The  decay  of  the  I-V  characteristics  from  their 
initial  to  final  state  was  easily  recognizable.  The  decay  corresponded 
to  a period  less  than  a second,  which  is  consistent  with  the  thermal 
time  constant  measurements  done  at  this  temperature. 

In  the  measurements  described  above,  the  decay  was  caused  by  the 
thermal  emission  of  holes  from  the  gallium  centers  to  the  valence  band. 
However,  should  the  surface  of  the  device  be  illuminated  by  infrared 
radiation  in  the  8-14  micrometer  wavelength  range;  and  of  sufficient 
intensity  to  optically  ionize  all  the  gallium  centers  in  the  depletion 
region,  then  the  magnitude  of  the  decay  at  24°K  would  be  identical  to 
this  shown  in  Figure  27 . The  decay  in  the  later  case  would  be  completed 
in  a shorter  period,  which  corresponds  to  around  four  times  the  optical 
time  constant.  In  other  words,  one  can  interpret  the  dark  and  the  dashed 
lines  in  Figure  27  as  the  I-V  characteristics  of  an  1RFET  before  and  after 
Illumination  respectively.  Consequently,  we  observe  that  the  current 
decreases  in  the  IRFET  after  illumination.  This  result  is  consistent  with 
the  prediction  made  in  the  IRFET  mathematical  model  (Appendix  Equations  24 
and  25). 


GALLIUM -BORON  DOPED  MOSFET 
Vbq=  0-0 
TEMP  * 24°  K 


■DARK  UNES- GALLIUM  NEUTRAL 


^DASHED  UNES -GALLIUM  NEGATIVE 
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8 

Vos  (VOLTS) 

FIG.  27  The  I-V  characteristic  of  an  IRFET 

before  and  after  illumination. 
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The  photocurrent  Al  Is  the  difference  In  current  at  two  corresponding 
1/5 

points  on  the  dark,  and  dashed  lines  of  Figure  27,  at  the  same  electric 
bias  conditions.  In  the  linear  region  of  the  I-V  characteristics  of 
Figure  27  , Alog  did  not  appear  to  be  fully  independent  of  the  gate 
to  source  voltage,  as  predicted  by  the  IRFET  mathematical  model,  (Appendix 
Equation  24).  However,  no  detailed  investigation  of  this  dependence 
in  the  linear  region  was  attempted,  therefore  one  would  not  attach 
great  importance  to  our  above  mentioned  observation. 

The  MOSFET  current-voltage  characteristics  in  the  saturation  region, 
(Appendix  A,  Equation  22)  was  evaluated  at  several  points  on  the  I-V 
characteristics,  at  various  back  gate  bias  potentials,  in  order  to 
determine  an  average  surface  electron  mobility  that  could  be  used  in 

2 

the  device  equations  to  model  the  IRFET.  Mobility  values  of  257.0  cin  /y. sec 

at  29.3°K,  286.4  cm^/V.sec  at  24°K  and  286.6  cm^/V.sec  at  23.09°K  were 

obtained  at  VGg  - VT  = 35  volts.  These  electron  mobility  values  are 

very  low  for  the  (100)  orientation  at  this  temperature  range.  An  electron 

2 

mobility  value  above  1200  cm  /V.sec  was  expected,  since  the  lattice 
scattering  has  been  greatly  reduced  at  this  temperature  of  operation  [ 29] . 
The  mobility  values  obtained,  also  show  very  little  sensitivity  to 
temperature  -1.66Z  per  °K.  These  observations  suggest  that  the  impurity 
concentration  is  very  high  near  the  surface.  If  we  start  with  the 
assumption  that  our  wafers  are  uniformally  doped  with  boron;  then  we 
have  to  conclude  that  the  gallium  penetration  was  not  5 micron  deep 

below  the  surface  as  expected.  Other  measurements  done  on  our  devices  were 
found  to  confirm  this  observation,  like  the  measurement  of  the  junction 
characteristics  obtained  after  gallium  diffusion  in  P-type  wafers  and  the 
dependence  of  AV.J.  on  the  back  gate  bias  potential. 


FIELD  EFFECT  MOBILITY,^  •(2/cy(L/W)I^(\^s-Vrf 
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This  conclusion  contradicts  some  of  the  results  obtained  by  others 
( 10  ] on  the  impurity  redistribution  in  silicon  due  to  thermal  oxidation. 

It  also  makes  it  clear  that  inspite  of  the  30  minutes  predeposition 
period  at  1200°C  and  the  60  minutes  drive  in  diffusion  at  1100°c  in  wet 
t^ie  Rallium  profile  is  far  from  being  uniform  in  our  devices.  More 
experimental  work  needs  to  be  done  on  the  gallium  diffusion  process,  in 
order  to  obtain  a lower  surface  concentration  and  better  surface  parameters. 

The  square  root  of  the  source  to  drain  current,  /log,  in  the  saturation 
region  was  plotted  versus  the  gate  to  source  potential,  VGS,  Figure  29  , 
with  the  back  gate  bias  potential  as  a parameter.  The  slopes  of  the 
resulting  lines  were  used  to  determine  the  previously  quoted  mobility 
values  as  follows: 

IDS-  U CQ  M ( VGS.-_VT)  2 

L 2 (Eqn,  (22),  Appendix) 

Therefore 


Therefore  the  slope  in  Figure  29  corresponds  to 


V^DS 
vGS-  VT 


- 2.3  x 10"-*  /Amp/volt. 


(86) 


(87) 


2 J_  2 

. . p » Cc  W/L  (slope  ) (88) 

where  W/L  is  the  channel  width  to  length  ratio  » 8.17  CQ  - 5.06  x lO^pf  = 
gate  capacitance  per  unit  area  based  on  7000A°  oxide  thickness. 

Therefore 

W * 257.06  Volt/cm^. second.  (gq) 


This  mobility  value  was  verified  by  several  other  measurement  sets 
at  29.4°K. 


- - mr 
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We  calculated  the  electron  mobility  from  equation  (88)  as  the 
temperature  went  down.  The  mobility  peaked  around  100°K  then  dropped 
slightly  in  value,  then  increased  to  a maximum  of  around  286  cm'Vvsec 
at  23°K.  The  results  of  these  calculations  are  plotted  on  Figure  28. 

F.  Fang  and  A.  Fowler  have  also  measured  |^2‘)J  the  effective  mobility  on 
the  (100)  orientation.  They  kept  however,  the  source  to  drain  voltages 
low  to  avoid  the  hot  electron  effect  and  the  non-uniformity  of  the 
surface  field  across  the  sample  | . They  reported  that  the  mobility 
measured  at  varying  electric  field  at  low  temperatures,  has  increased 
initially  to  a 1300  cm?/Vsec  peak  at  2 x 10*J  V/cm  electric  field  then 
dropped  rapidly  as  the  electric  field  was  increased  further.  The 
samples  used  in  these  measurements  were  characterized  by  high  substrate 
resistance  or  very  low  doping  (129  ohm-cm) . 

The  mobility  has  increased  consistently  as  the  temperature  went 
down.  The  mobility  peak  (1300  cm^/Vsec)  was  obtained  at  around  30C>K, 
thereafter  much  less  temperature  sensitivity  was  reported.  This  last 
set  of  measurements  was  carried  out  on  a 1.88  ohm-cm  sample,  (100) 
orientation  and  when  the  gate  voltage  was  15  volt  above  the  threshold. 

These  results  are  consistent  with  the  general  trend  that  we  have 
observed.  The  deviation  is  in  the  order  of  magnitude  of  the  surface 
mobility,  as  a consequence  of  the  gallium  accumulation  at  the  surface. 

We  found  also  that  there  is  no  single  power  law  exp  that  can  explain 

the  mobility-temperature  dependence  over  the  entire  observation  range. 

The  temperature  dependence  of  the  threshold  voltage  on  device  A 
was  measured.  On  Figure  28  we  plotted  the  shift  in  the  threshold 
voltage  from  its  value  at  room  temperature  versus  1000/T.  The  shift  in 


■■■I 


GALLIUM  DOPED  MOSFET  - SATURATION  REGION 
DEVICE  A Nj  = 2.25  x I0,6cm2 
TEMP  * 29.3°K 


GATE  VOLTAGE  V„e(  VOLTS) 

6S 

FIG.  29  A plot  of  drain  current  versus  gate 
voltage  in  the  saturation  region 
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the  threshold  voltage  as  the  temperature  is  lowered  was  also  observed  by 
C.  Rogers  {joj  , F.  Fang  and  A.  Fowler  [j?9j  . The  threshold  voltage 
is  given  by  equation  (6)  in  the  Appendix  ^27j  . From  this  equation 

we  should  conclude  that  the  primary  dependence  of  the  threshold  voltage 
on  temperature  results  from  Fermi  level  dependence  on  temperature.  The 
Fermi  potential  appears  as  a part  of  the  surface  potential  <j>s , and  the 
work  function  $ms , in  the  above  mentioned  equation. 

If  <J>S  is  assumed  to  be  twice  the  Fermi-potential,  and  if  the 
dependence  of  $ms  on  the  Fermi  potential  is  added,  then  we  should 
expect  that  the  shift  in  the  threshold  voltage  with  temperature  should 
not  exceed  three  times  that  in  the  Fermi  potential  or  in  other  words 
it  should  not  be  more  than  one  volt  at  the  most.  However,  the  observed 
shift  in  the  threshold  voltage  is  much  higher  than  that  (above  20v 
between  300°K  and  30°K) . This  has  also  been  observed  by  many  researchers, 
for  example  ^29j  , and  the  shift  was  attributed  primarily  to  the 

surface  charge  in  the  oxide  close  to  the  Si-SiOp  interface.  Therefore 
we  will  proceed  to  determine  the  change  in  the  surface  charge  w.r.t.  a 
typical  value  of  say  Qss/q  - 9 x 10^°/cm2,  that  could  have  caused  the 
observed  shift  in  the  threshold  voltage.  The  procedure  for  doing  this 
is  simply  to  determine  the  threshold  voltage  at  the  above  mentioned 
typical  value  for  surface  charge,  then  to  match  this  with  the  measured 
threshold  voltage  and  hence  determine  the  deviation  of  the  actual 
surface  charge  from  this  typical  value.  The  results  of  our  calculations 
are  plotted  on  Figure  28. 

Another  low  temperature  observation  was  the  hysteresis  appearing  in 
the  drain  characteristics  on  the  curve  tracer  near  30°K  and  low  current 
values.  The  observed  hystero-iis  virtually  disappeared  when  the  gate  voltage 
was  doubled.  The  hysteresis  effect  was  attributed  to  trapping  [30)  mechanisms. 
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We  will  not  attempt  to  analyze  tills  observation  in  any  detail  here. 

When  the  IRFET  was  triggered  by  a square  wave  on  the  gate,  at 

low  temperatures,  the  rise  time  of  the  current  between  tin-  source 
and  drain  was  short  and  temperature  independent.  The  current  reached 
the  value  I ^ almost  instantaneousl yt  with  the  rising  edge  of  the 
gate  voltage.  When  the  TRFF.T  was  turned  on  with  the  positive  gate 
voltage,  the  source  to  drain  current  decayed  with  the  thermal  time 
constant, to  it's  final  value  I A voltage  proportional  to  the  source 

to  drain  current  was  displayed  on  the  oscilloscope.  This  voltage 
was  dropped  across  a 96.8  ohm  resistor  at  the  input  of  the  oscilloscope. 

The  currents  I ^ and  Al  _ were  measured  on  Device  A,  N = 2.25  x 

Uo  Uj  1 

16  -3 

10  cm  , at  various  gate  to  source  voltages.  The  source  to  drain 
voltage  was  kept  constant,  at  V = 10  V.  The  measurements  were  repeated 
at  different  temperatures  and  at  different  source  to  drain  voltages. 

The  IRFET  photocurrent,  -Al^ , is  shown  plotted  in  Figure  30,  as 
a funtion  of  V in  the  saturation  region,  at  V = 10  V,  V =0,  and 

Vjo  1)o  Ijv» 

at  29.4°K.  The  plotted  line  corresponds  to  the  function  predicted  by 
equation  (25)  in  the  attached  appendix. 

Calculations  based  on  the  measured  photocurrent,  -AlpC , and  the 
measured  source  to  drain  current,  I *1  , at  29.4°K  and  V = 10  volt, 
enabled  us  to  plot  the  dependence  of  the  initial  source  to  drain  current 
on  Vgg , Figure  31.  This  graph  shows  the  dependence  of  the  initial  and 
final  source  to  drain  currents  on  V , in  the  saturation  region.  The 
two  lines  in  Figure  31,  upper  and  lower,  represent  the  two  charge  states 
of  the  device,  gallium  neutral  or  negative  respectively.  AVT  was 
determined  from  Figure  31,  it  corresponded  to  0.9  volt. 
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Figure  32,  shows  the  threshold  voltage  dependence  on  the 

/I  + V , where  V is  the  back  gate  bias  potential.  The  threshold 
BO  BO 

voltage  values  were  determined  at  the  points  of  intersection  of  Mtc 

lines  with  the  V axis  on  Figure  31  . 

Lib  GS 

Figure  33  , shows  the  change  in  the  device  threshold  voltage  bet- 
ween the  two  charge  states  as  a function  of  the  back  gate  bias  potential 
Equation  26,  Appendix,  was  used  to  calculate  the  change  in  the 
threshold  voltage  from  the  measured  data,  at  V = 0,  as  follows: 


V - V 
( GS  T) 


= (35/2)  6.5 


- 1.2  volts 


The  total  substrate  doping  for  device  A was  determined  from  equation  (5) 

and  Figure  (19),  N , * N,  - 2.25  x 1016cm~3.  The  Ga  & B concentrations 
sub  I 

will  be  calculated  from  equation  18,  Appendix,  as  follows: 

d (v.)  r i / \ I 


4V*C  UBG  • °> 


ion  io,  Appendix,  as 

(V  f-i-l  "OA 
'i*T  L2  ( Ni  I . 


Therefore; 


1-2  ■ <12>  <i>  ohrro16 


Therefore ; 


Therefore ; 


0.45  x lO^cm  3 


na  “ Ni  - nga 


2.25  x 1016  - 0.45  x 1016 


1.8  x 1016  cm"3 


The  theoretical  value  of  the  threshold  voltage  of  the  device  was 
determined  from  the  well  known  MOSFET  equations,  equations  (6)  through 
(11)  in  the  Appendix;  the  following  discussion  illustrates  the  methods 
we  used  in  this  calculation 
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FIG.  32  The  dependence  of 
threshold  voltage  on  back-gate 
bias  at  low  temperature. 
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The  effective  density  of  states  in  the  valence  and  conduction 


bands,  N and  N were  calculated  at  29.4  K: 

V Vi 


2n  m K.T  ' 
NV  “ NC“[  ^2  * 


Therefore , 


19  29  4 

• Nc  ■ Nv  " 2’5  x 10 


7.669  x 1017/cm3 


under  the  assumption  that  the  electron  and  hole  masses  are  equal  to 
the  mass  of  the  free  electron: 

me  ’ "*h  “ mo  (97) 

The  silicon  bandgap  varies  with  temperature  [ 31  ],  according  to 

the  following  equation 

oti 

E (T)  - E (0)  - T + 6 (98) 

8 8 

Eg  (0)  was  measured  [ 32  ] , and  found  equal  to  1.16  eV  at  0°K. 
The  coefficients  a and  0 were  also  determined  [ 31  ] for  silicon,  and 
were  found  to  have  the  following  values 


a - 7 .0  2 x 10 


(100) 


0 - 1108 


(101) 


The  difference  in  thermal  expansivity,  between  the  header  and 
the  device  should  cause  pressure  that  may  alter  Eg.  However,  this  was 
neglected  in  our  calculations. 


E (29.4  K)  - 1.16  - 

8 


7.02  x 10~A(29.4)2 
29.4  + 1108 


= 1.16  eV 


(102) 


The  intrinsic  concentration  of  electrons  in  the  conduction  band  at  29.4  K 
was  determined  from  the  ‘ollowing  equation: 


"l  “ /NcN;  exp  1 kt”1 


(103) 
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4 


Therefore , 


V2,-4°K)  - (7.669  * 1017)  oxp  f-  oioOZSZS1 


4X10  182/cm3 


(104) 


as  compared  to  1.4  x 103^/cm3  at  300°K 


where 


29.4 

— (at  29.4°K)  - — (at  300°K)  "l00 


— “ °inn~  29-A  = 0.002528  (105) 

qjUU 

KT  P 

The  Fermi  level  (U„|  * — In  ) , at  the  temperature  of  operation 
r q n . 


lies  between  the  gallium  and  boron  levels.  Therefore: 


4>F  = ^ 


-.  58eV 


(106) 


At  this  low  temperature  range,  around  29.4°K,  all  the  boron  impurity 
atoms  were  still  ionized.  Therefore,  the  carrier  concentration  during 


the  read  in  period  is  given  by: 

P - N ~ N + N * 2.25  x 1016cm'3. 
1 A (jA 


(107) 


The  nobility  is  known  to  decrease  with  the  high  electric  field, 

3 

where  the  high  electric  field  operation  is  defined  as  E > 10  V/cra. 

This  is  known  as  the  "hot  carrier  effect".  The  mobility  dependence 
on  temperature  and  orientation  is  also  complex.  Therefore  we  have  used 
the  measured  mobility  value  at  this  temperature  in  the  model  calculations. 
U e (29.4°K)  - 257.06  cm2/V.sec.  (108) 

The  work  function  between  the  semiconductor  body  and  the  pate  field 
plate  was  calculated  for  the  A1  - S102— Si  system,  as  follows: 


>m  “ " *mo  " (*so  + 2q  + *F> 


(109) 


3.2  - (3.25  + + 0.58) 


-1.21  volt 


(110) 


p 
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where  fc  and  $ are  the  barrier  lirtiftlir « between  the  conduction  hand  of 
IDO  so 

the  SiO^  and  both  the  fermi  level  of  the  gate  metal, and  the  silicon 
conduction  band  respectively. 

The  surface  charge  Qsg  at  the  Si  - SiO^  interface  is  strongly  process 

dependent  and  may  vary  with  temperature.  We  have  assumed,  however,  the 

following  typical  value: 

Q 

— - —9  x 1010/cm2  (111) 

q 

The  gate  oxide  capacitance  was  calculated  from  the  knowledge 
of  the  oxide  thickness  (7000A°). 

Co  - Vox  - 5.057  X 10-9  F/cm2  (112) 

^OX 

The  ref  ore 


-Qss  " (9  x IQ10) (1.6  X 1019)  = -2.84  volt 

(T-  5.057  x 10-9 


(113) 


The  flat  band  voltage  can  be  determined  as  follows: 

Qss 

vr„  * (i  - tt~  - -1.21  - 2.84  « -4.05  volts  (114) 

r O TUS  LO 

From  equation  (10)  in  the  Appendix,  we  determine  (A), 

A - VpB  + 2<{>f  « -4.05  + 1.16  - 2.89  volt  (Hr>> 

Similarily,  (B)  can  determined  from  equation  (11)  in  the  Appendix: 


B ■ 'S,i  Co"  <2»f  - V> 

Co  

- ^X  12  X 8.85  K IQ-1**  x 1.6  x 10~19  x 1.16 


5.057  x 10 


(116) 


1.24  x 10  2 (at  zero  back  gate  bias)  (117) 

The  initial  value  of  the  threshold  voltage,  follows  directly 
from  equation  (12)  in  the  Appendix. 

V*  - A + B (118) 

- -2.89  + 1.24  x 10_7(1.8  x 1016)*5 


* 13.7  volt 


(119) 


The  final  value  of  the  threshold,  follows  from  equation  (13)  in  the 
Appendix: 
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VT  " A + B /nA  + Nga  (120) 
= -2.89  + 1.24  x 10" 7 1 2 . 25  x 101{' 

=»  15.7  volt  (121) 
now  if  we  match  V*  as  calculated  from  equation  (119)  to  that  measured 
we  can  determine  the  actual  surface  charge. 


-Qss 

Co 


V 


1 

T measured 


+ 


9 x 1010 

q 


(122) 


The  initial  and  final  threshold  voltages  as  determined  at  29.4°K  by 
extrapolation  on  Figure  31  were  4.7  and  5.8  volts  respectively. 
Therefore 

"aT  ” 4‘7  " ( 13-7  + 2-84  ) (123) 

and  hence  therefore 

^ - 3.74  x 101 1/cm2  (124) 

9 

at  29.4°K 


Similarly,  the  actual  surface  charge  was  determined  at  other 
temperatures  and  the  result  was  plotted  on  Figure  28. 

The  threshold  voltage  that  was  originally  negative  at  room 
temperature,  changed  to  be  a positive  quantity  at  low  temperatures, 
i.e.  the  depletion  device  changed  to  be  an  enhancement  device  at 
low  temperatures.  No  significant  changes  in  the  threshold  voltage 
were  noticed  in  the  temperature  range  between  21°K  and  35°K. 

On  Figure  33  we  notice  that  the  magnitude  of  AVT  dropped  as 
the  back  gate  bias  potential  was  increased.  This  is  contrary  to  the 
behavior  that  our  model  suggests.  However,  the  model  is  only  valid 
for  uniformly  doped  devices.  The  gallium  accululation  near  the  surface 
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results  In  a drop  in  the  doping  ratio  of  the  depletion  region  N_./N., 

1>A  A 

as  the  size  of  the  depletion  region  Increases  with  the  back  gate  bias. 

This  fact  ultimately  leads  to  the  reported  decrease  in  AVT  as  the  back 

gate  bias  potential  is  Increased. 

The  change  in  the  threshold  voltage  AV.^,  was  calculated  from 

equation  (25)  of  the  Appendix  at  V„c  - 40,  V*  + AV/2  - 16.5V  (Figure  30), 

05  T T 

W/L  » 8.17,  u - 257.06  cm2/V.sec,  Co  ■ 5.057  x 10  °F/cm2  and 
A1DS  - 0.43  mA. 


'AIDS  * " (AV  ( VGS  - VT  - AVT/2 


(125) 


Therefore 

AVt  - 1.2  volt  (126) 

This  compares  very  closely  to  the  value  of  AV^  determined  by 
extrapolation  of  the  experimental  data  on  Figure  31. 


C.  Optical  Mt- ns i iremon l s 

When  tlie  temperature  Is  around  16°K,  the  time  constant  of  thermal 
ionization  of  the  gallium  centers  is  several  seconds  long.  The  optical 

emission  rate  at  this  temperature  range  is  much  higher  than  the  thermal 

. . o t 

emission  rate  c >>  e . 

P P 

Where  the  optical  emission  rate  is  defined  as 

e°  - o°$  (127) 

P 

o°  is  the  photo ioni zat ion  cross  section,  which  is  a measure  of  the 
probability  of  carrier  emission  when  photons  of  a certain  frequency 
are  incident  on  the  device,  and  $ is  the  incident  photon  flux  in  number 
per  unit  area  per  second. 

The  target  of  the  optical  measurements,  is  to  verify  that  the 
Ga-doped  MOSFET  can  operate  as  an  infrared  detector  for  radiation  in 
the  far  infrared  range  (8  - 14  ura) , as  well  as  to  measure  the 
photoionization  cross  section  of  the  gallium  centers  at  various 
frequencies  in  the  far  infrared  radiation  range. 

To  separate  the  optical  response  from  the  thermal  response,  the 
photoionization  cross  section  has  to  be  measured  at  temperatures  around 
or  below  21°K.  At  this  temperature  the  thermal  decay  is  negligible. 

Our  thermal  decay  measurement  reveals  that,  the  activation  energy 
(E  - E^)  of  the  Ga-Center  has  strong  eLoctric  field  dependence. 

Ega  - Ey  can  be  varied  from  0.071  eV  at  a field  strength  of  10?V/cm 
to  0.027  ev  at  3 x 1 0 '*  V/cm  which  is  quite  different  than  0.065  eV 
given  by  ref.  (lft).  If  the  photoionization  cross  section  has  the 
same  field  dependence,  it  is  conceivable  that  the  spectral  response 
of  the  Ga-doped  MOSFET  can  be  extended  well  into  the  infrared  by 
applying  a strong  electric  field.  In  other  words,  a tunable  infrared 
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detector  can  be  realized. 

In  our  experiments  - so  far  - we  have  aimed  at  a gallium 
concentration  lower  or  much  lower  than  the  boron  conoenLrnt ion 
in  the  substrate.  However,  better  optical  detection  and  device 
performance  should  be  expected,  when  the  gallium  concentration  is 
much  higher  than  the  boron  concentration  in  the  substrate.  This 
will  obviously  increase  the  photo  current,  Al  ,,,  appreciably. 

Although  the  IRFET  can  be  operated  in  both  the  linear  and  satu- 
ration region  the  responsivity  is  obviously  higher  in  the  saturation 
region.  To  illustrate  the  integration  function,  which  is  a property 
of  the  IRFET,  the  photocurrent  time  dependence  in  the  saturation  region 
is  described  in  the  following  equation,  which  compares  to  equation 

(25)  in  the  Appendix  when  AV  /2  <<  V - V 

T GS  T 


AIDS(t)  = ~i^L  (AVt)  (V(;s-V,j)  (l-exp(-t/T>) 


(128) 


If  the  time  period  t <<  x then  equation  (128)  can  be  simplified 
to  the  following 

4IDS(t)’,'lIDS  7 <>29) 

Therefore,  the  IRFET  integrates  the  incident  photon  radiation  in 
the  time  period  t <<  t.  The  total  input  energy  is  therefore  expressed 
as  follows  (7j 

E.  * <t.  t.  tiw.  A (130) 

in 


where 


^ - , where  h is  the  Plank  constant  = b.625  x 10-3,‘  Joule,  second. 

ZTI 

w « 2tiv  where  v is  the  frequency  of  the  incident  radiation. 

A is  the  active  area  of  the  device  ■ WL 


The  output  power  at  the  end  of  the  integration  period  t (t<<x),  is 


/(> 


Dj Per  for mane o 

In  this  section  we  wiil  attempt  to  calculate  a typical  value  for 

the  responsivity  of  a Ga-B  doped  IRFliT.  The  device  will  he  assumed 

to  have  the  following  specifications 

W/L  = 8,  A = 2.5  x 10”  W,  X = 7000  A° 

o 

The  input  energy  can  be  calculated  from  equation  (I'JO) 


E.  =$ . t .■fin). A , 
in 

The  photon  energy  at  II  n meter  is  given  by 
4*  . . >•  6.25  x 10- x 3 x I010 

1.6  x 10“*^  1 x I0-1' 


if 


♦ = 1016  11/ cm'  sec  and  t = 500  mi  1 1 i sec 

Then  the  input  energy  can  be  calculated  as  follows 
E 


in 


1016  x 0.5  x -.107  x 1.6  x I0-10  x 2.5  x 10" 3 

= 2.14  x I0-7  Jouies 
The  output  power  = |Alpgl  x V^g 

= 0.55  x 10  = 5.5  miili  Watts 

Therefore  the  responsivity  can  be  determined  as  follows 


_ . Output  Power  5.5  miili  Watts 

Responsivity  = J - . — , 

Input  Energy  0.214  micro  Joule 


(134) 


(135) 


(136) 


(137) 


VI.  CONCLUSIONS 


u 


/ / 


A process  for  fabricating  the  gallium-boron  doped  MOSFET  has 
been  developed.  The  gallium  was  diffused  in  the  wafer  through  a 
silicon  dioxide  layer  (33). 

The  ionization  energy  of  the  gallium  center  was  measured  at  low 
temperature  in  the  depletion  region  of  a MOSFET  behind  a strongly 
inverted  channel.  The  thermal  activation  energy  of  the  gallium  center 
was  measured  and  found  less  than  the  previously  reported  ionization 
energy  of  gallium  [ 16, 19,2()|".  The  reduction  in  the  ionization  energy 
was  attributed  to  the  barrier  lowering  as  a result  of  the  MOSFET's  surface 
potential.  The  Poole-Frenkel  model  [9]  was  used  in  the  data  analysis. 

It  was  found  that  the  barrier  lowering  calculated  from  this  model  is 
comparable  to  the  difference  between  the  previously  reported  gallium 
ionization  energy  (19]  and  the  measured  thermal  activation  energy. 

A solution  of  Poisson's  equation  in  the  depletion  region  of 
the  MOSFET  suggests  that  part  of  the  gallium  population  in  the  depletion 
region  may  be  permanently  ionized  at  any  low  temperature  by  the  electric 
field. 

The  mobility  and  threshold  voltage  temperature  dependance  that 
were  observed,  were  consistent  with  the  low  temperature  observations 
of  others  [29,30].  The  observed  shift  in  the  threshold  voltage  was 
attributed  to  changes  in  the  surface  charge. 

The  1 Rl'ET  model , previously  suggested  [2]  was  found  to  describe 
the  C.a-doped  IKFET  operation  with  accuracy  sufficient  for  most 
engineering  applications. 
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APPENDIX 

W.  C.  I’arkor  ami  I.  Forbes  1 have  presented  a model  of  t lie  I KI'ET 
operation.  The  gallium-boron  double  doped  MOSFKT  responce  data  was  found 
in  quite  good  agreement  with  that  predicted  by  theory.  The  model  equations 
are  summarized  here  for  easy  reference.  The  equations  were  also  modified 
to  fit  the  operation  of  the  double  doped  gallium-boron  N-channel  MOSFET. 

To  explain  the  phenomena  of  the  modulation  of  the  channel  conduc- 
tance by  the  infrared  radiation  rt  low  temperatures  mathematically,  we 
will  start  by  stating  the  charge  neutrality  condition, 

^G  + QSS  = " ^QCH  + ^Ga  + V ^ 

where  Q is  the  charge  on  the  gate,  Q is  the  fixed  positive  surface 
p."  state  charge  at  the  silicon-silicon  dioxide  interface,  is  the  charge 

contributed  by  carriers  in  the  channel,  Q is  the  charge  of  ionized 
gallium  atoms  in  the  oxide  or  near  the  silicon  surface,  and  is  the 
charge  in  the  surface  depletion  region. 

When  the  gallium  centers  in  the  depletion  region  are  photoi onized , 
the  net  space  charge,  Q , becomes  more  negative,  and  the  MOSFET  has  to 

D 

lose  electrons  from  the  channel  to  reduce  the  channel  charge,  Q , and 

‘UK 

maintain  neutrality.  This  results  in  a corresponding  decrease  in  the 
magnitude  of  the  drain  current.  All  other  terms  in  equation  (1)  are  con- 
stants . 

1 . The  Time  Constants  of  the  Device 

The  thermal  time  constant  will  be  referred  to  as  T , and  will  be 

defined  as  the  reciprocal  of  the  thermal  emission  rate,  cC.  Sim- 

P 

ilarly  the  optical  time  constant  will  be  referred  to  as  r , and  will 

o 

be  defined  as  the  reciprocal  of  the  optical  emission  rate,  e°. 

P 


t 


I 


n 


(2),  O) 


2 . The  Effective  Doping  of  tlif  I'-'l'ype  Substrate 

Nj  is  the  effective  doping  of  the  I’-type  substrate,  is  the  boron 

concentration,  and  N„  is  the  pallium  concentration. 

(»a 

(a)  During  the  reset  operation  with  the  accumulation  voltage  applied: 
Since  the  gallium  centers  are  neutral,  then 


N, 


N. 


(M 


I A 

(b)  The  condition  during  the  read-in  period  or  the  device  effective 
doping  after  the  photoionization: 

Since  the  gallium  centers  are  negative,  then 

N.  ■ "a  " "c.  <5) 

3 . The  Threshold  Voltage 

For  an  n-channel  MOSFF.T,  the  threshold  voltage  is  given  by  ^27) 

qr 

V = V + 2<i>  + — (h) 

T FB  VF  c ^ 

O 

where  V is  the  gate  voltage  required  to  establish  the  flat  band 
r n 

condition,  <t>p  is  the  Fermi  potential,  is  the  gate  oxide  capaci- 

tance/unit area,  and  Q is  the  depletion  layer  charge  generated  hv 

D 

the  ionized  impurity  atoms. 


% ' - VBR> 


(7) 


Wliere  K is  the  rel.ative  dieiectric  constant  of  silicon,  e is  the 
Si  o 

electric  permitivity  of  free  space,  q is  the  electronic  charge,  and 
VBG  18  the  back  gate  bias  potential:  a reverse  bias  on  the  substrate 

that  increases  the  electric  field  in  the  depletion  region  and  ionizes 
additional  impurity  centers. 


- 


*(  I 


substituting  equation  (7)  into  equation  (6),  we  obtain 

V,.  - VrlJ  + 2$,  + 


FB 


^SlV'^F  - V *r- 
■t  ■ — — c /NI 


(8) 


Equation  (8)  can  also  be  expressed  as  follows: 


V - A + B 


(9) 


where 


A - Vre  + 2V 


(10) 


and 


^c7c  nl2^  - v'  ] 


Si  o"  7 BO 


(11) 


A and  B will  be  assumed  constants  for  a device  at  constant  temperaUre 
and  constant  back,  gate  bias  operating  conditions.  The  only  term 


that  may  vary  is  since  it  is  dependent  on  doping.  However,  for 


V > 2<,>  our  assumption  will  still  be  accurate. 

Dlf  I 


(a)  The  threshold  voltage  during  the  reset  operation  with  the  ac- 
cumulation voltage  applied  (initial  condition); 

(12) 


v|  * A + B /Na 


as  obtained  from  substituting  Equation  (4)  into  Equation  (9). 

(b)  The  threshold  voltage  after  tie  photo ionizat ion  (final  condi- 
tion) ; 

by  substituting  Equation  (5)  into  Equation  19)  we  gel ; 


V*  *=  A + H «4T  + N 
1 A (jrt 


(ID 


Hence,  the  increase  in  the  threshold  voltage  is 


AV 


i ' < - \ > 


(14) 


~ V *+ 


, / " «Cn 

B H 
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(15a) 


for  tho  a iso  that  N,,  r-  N ; wi* 
(..I  A 

AV  iwfi7(  I + (N  . /2NJ  - I)  - 
1 A (..i  A 


( I M. ) 


(c)  The  variation  of  the  threshold  voltage  ns  a function  of  the 
back  gate  bias  potential; 

differentiating  equation  (9),  under  the  assumption  that  24' j- 


d(VT) 

d[/TTT\ 


Ivbg|] 


b 


(K>a) 


B = — [d(V  )/d(/l  + V )] 

^ T BC 


(lhb) 


Equation  (15b)  for  the  case  that  N <*  N^,  can  be  expressed  in 
terms  of  the  effective  doping  N as  shown  in  equation  (17) 


2 ^ 


Substituting  equation  (16b)  into  equation  (17),  and  assuming 

N„  <<  N,  we  obLnin 
Ga  A 

d(V  ) N 

av  - I I _££ 

T di/TT*-|  2 "i 

In  other  words,  the  gallium  concent t at  loo  in  the  device  can  he 
determined  from  the  measurement  of  AV^,  at  a certain  back  gate  bias 
potential  and  from  the  measurement  of  the  V.^,  depcmlanco  on  the  Back 
gate  bias  potential. 


tY 


Tin*  Time  Dopeinlnnce  of  the  Ionization  Process : 


The  number  ol  Ca 


centers  in  the  negative  charge  state  N is  given  by: 

-t/T 

NT(t)  * NCa(l  - e °]  (19) 

under  the  assumption  that 


r '•>  T at  the  temperature  of  operation, 
t o 

Where  N^.(t)  is  the  total  number  of  the  ionized  gallium  centers  at 
t ime  t . 

The  Time  Dependance  of  the  Threshold  Voltage 
From  equation  (15-a)  and  (19),  we  obtain 


= V*  + AVT(t) 


vi  +h/h7[A  + rr3-  [i  - 


-t/T 


11 


(20) 


where  V^(t)  is  the  threshold  voltage  as  a function  of  the  illumina- 
tion time. 

MOSFF.T  Current-Voltage  Equation  in  the  Linear  Region 

The  linear  region  is  defined  in  the  well  known  [27]  MGSFKT  model  as 

the  region  of  the  1-V  Characteristics  where 


V - V > V 
GS  T DS 


Therefore  the  source  to  drain  current  I in  the  linear  region  is 


DS 


VC  W 

2L  t2VDS(VGS"  VTJ  ~ VDf^ 


(21) 


W 

where  p is  the  carrier  mobility  at  the  surface,  — is  the  channel 
width  to  length  ratio,  V^s  is  the  drain  to  source  voltage,  and 
is  the  applied  gate  voltage. 


HOSFET  Current  Voltage  Equation  in  the  Saturation  Region 

For  V > V - V , the  MOSFET  is  operating  in  the  saturation 
DS  Gs  T 

In  this  case  (27) 


region. 


H(> 


nC  w 

I = — — (V  - V )“ 
l)S  2L  V CIS  V 


C in  equations  (21)  and  (22)  is  the  oxide  capacitance  per  unit  area, 
o 

C.  can  be  determined  from  the  parallel  plate  capacitance  equation, 
o 

C - k r /t  (23) 

o ox  o ox 

where  k is  the  dielectric  constant  of  the  oxide  and  t is  the 
ox  ox 

oxide  and  t is  the  oxide  thickness  under  the  gate, 
ox 

8.  The  Photocurrent 

The  photocurrent  Al^,.  is  the  change  in  the  source  to  drain  current 
as  a result  of  the  photoionization.  Al  can  bo  expressed  in 
terms  of  the  change  in  the  threshold  voltage  AV  for  both  the  linear 
and  saturation  regions. 

Foi  the  linear  region  from  equation  (21),  for  a given  VDS 


-pC  W 
o 


(AV  ) (V  ) 
v T DS' 


For  the  saturation  region  from  equation  (22) 
-pC  W 

ALDS  " L (Avt>  (Vgs  - VT  - AVt/2) 


The  negative  sign  in  the  two  above  equations  indicates  that  the  source 
to  drain  current  decreases  after  illumination.  Dividing  equation 
(25)  by  equation  (22),  we  have: 

.1  , (2  AV  ) 

(2b) 

DS  V ;s  - VT 

in  the  saturation  region. 

9 • The  Current-Time  Dependance 

From  equations  (20),  (24),  and  (25)  we  obtain  for  the  linear  region: 


Hi 


+IIC  w 
o 


l>S 


(vgs  - 


v.r  - b*4t  |/ 

A 


A + 


n p - exp  -t/i  i 

(•«*)  o 


(27) 


For  the  saturation  region: 


■ds"’ 


+uc  w 

o 

2L 


IV, 


CS 


- vi  - b/n~  [/i  + 

i A 


N[l-  exp 

u3 


-t 


TO 


] 


-HI  (28) 


Equations  (27)  and  (28)  express  the  source  to  drain  current  as  a 
function  of  the  illumination  time.  The  time  t = 0,  corresponds  to 
the  beginning  of  the  ionization  process. 


10.  The  time  dcpendance  of  the  threshoid  voltage  and  drain  current  when  the 
NCa  ~ NA  & ATbS  “ W 

If  NCa  <<:  N^,  equation  (20)  can  be  simplified  using  tiie  binomial 
approximation . 


t . BN„  . 

v!'  (t)  = vi;  + — — (1  - e_t ^°) 

2/N~ 

A 


(29) 


From  equation  (21)  and  equation  (24),  by  algebraic  manipulation 
we  obtain  the  following  equation 
Al. 


A V. 


US 


(V  - V - V ) 
1 j)s  GS  T US/2J 


(30) 


in  the  linear  region  ami  similarly  from  the  equations  (22)  and  (2b) 


wo  obtain 


88 


avt  • (vos  - V <'  - 1,1  + AW'ds  1 


(31) 


in  the  snturn L ion  region. 


II  A I |^j.  <v;  then  equation  (31)  c;m  l>o  simplified  i (>  i lu* 


following 


AI 


DS 


AV„.  = — — (V  - V ) 

i 2ius  cs  V 


(32) 


again  in  the  saturation  region. 


Equations  (27)  and  (28)  can  in  turn  be  simplified  when  N « N 

Ga  A ’ 

they  lead  to  equations  (33)  and  (34)  respectively. 


TDS(t> 


MCW  N 

~ VDS(VGS  ' Vx  “ B ~ d ~ exp  (~t/lo)) 

A 


V2> 


uc  w 

o 


2L  V'GS 


(V, 


BN 

" vt  “ ITiT  (1  " exp  c-t/ to) ) 2 

A 


(33) 

(34) 
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